U.S.  DEPARTMENT  OP  COMMERCE 
National  Technical  Information  Service 


AD-A028  512 


Statistical  Model 
of  Sonic  Boom 
Structural  Damage 

Booz-Allen  Applied  Research,  Inc. 


July  1976 


AO  A 0 2 8 5 l 2 


237103 


Itprt  Ik.  FM-H-TS-n 

STATISTICAL  MODEL  OF 

SONIC  BOOM  STRUCTURAL  DAMAGE 


Robert  L Hershey 
Thomas  H.  Higgins 


t 


July  1976 
Final  Report 


Document  is  available  to  the  public  through  the 
National  Technical  Information  Service 
Springfield,  Virginia  22161 


/DDC 


Prepared  for 

U.S.  DEPARTMENT  OF  TRANSPORTATION 

FEDERAL  AVIATION  ADMMISTRATION 
Sys.ems  Research  A Development  Service 
Washington,  D.C.  20590 


WHIOOUCtO  Bt 


NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

u s DffAtntiMT  of  cwetKC 
smiwcfiuc  VH  £2,41 


ir^f-  XWSS«*»»Wl3® 


1*  R*p«rf  N*.  2.  G*v»mm»nf  Acc#**i»n  N*. 

FAA-RD-76-87 

4.  Till*  «ii4  Subiitl. 

STATISTICAL  MODEL  OF  SONIC  BOOM 
STRUCTURAL  DAMAGE 

7.  A u*hot(t) 

Robert  L.  Hershey  and  Thomas  H.  Higgins 

9.  Performing  Organisation  Nome  end  Addrast 

Booz,  Alien  Applied  Research 
4733  Bethesda  Avenue 
Bethesda,  Maryland  20014 

32.  Sponsoring  Agency  Nome  end  Address 

Department  of  Transportation 
Federal  Aviation  Administration 
Systems  Research  and  Development  Service 
Washington.  D.C.  20590 

IS.  Svpplemenfery  Notes 


TECHNICAL  REPORT  STANDARD  TITLE  PaGF 

Fs.  Recipient* s Cetolog  No”  I 


5.  Report  Dote 

July  177 

A,  Performing  Organisation  Cod* 

I.  Performing  Orgenisetion  Reporr  No. 


10.  Work  Unit  No. 

11.  Confrect  or  Grant  No. 

DOT -FA  72  WA -2  823 

Type  of  Report  end  Period  Covered 

Fi  ial 

14.  Sponsoring  Agency  C©4<*> 

ATsT-y'* 


ia 


14.  Abtf'.ct 


The  probabilities  of  structural  damage  from  sonic  booms  were 
estimated  for  various  susceptible  structural  elements  using  a statistical 
modeling  technique.  The  breakage  probabilities  were  found  to  vary 
widely  with  the  specific  material  configuration,  but  to  consistently 
increase  with  increasing  nominal  overpressures.  The  ranges  of 
breakage  probabilities  at  a nominal  overpressure  of  1 psf  for  typical 
configurations  of  susceptible  materials  were  as  follows:  windows  — 

4 x 10"®  to  3 x 10“8,  plaster  — 3 x 10"^  to  5 x 10"?,  bric-a-brac  — 

1 x 10"®  to  1 x 10"8,  brick  walls — 6 x 10"?  to  1 x 10~9.  The  results  of 
the  modeling  tend  to  agree  well  with  sonic  boom  claims  experience. 


17.  IC.yV.r4:. 

Sonic  Boom 
Plaster 
Brick 
Glass 

IV.  $•« ytify  CImiiI.  ' *1 

Unclassified 

Fern  DOT  F 1700.7  ie.»> 


II.  Diitnkihtii  S'tttmMt 


Statistical  Simulation  Document  is  available  to  the  public  through  the 
Random  Loading  National  Technical  Information  Service, 

Bric-a-Brac  Springfield,  Virginio  22161 . 

Impact 


20.  foewiry  Citnil.  (•(  tk«»  p«go) 

Unclassified 


21.  N*.  .1  P I 22.  Pfij. 

U*~  \t.ocjlis 

FRKB  SUMECTTO  CHANGE 


The  contents  of  this  report  reflect  the  views  of  Booz,  Allen  Applied 
Research  which  is  responsible  for  the  facts  and  the  accuracy  of  the 
data  presented  herein.  The  contents  do  not  necessarily  reflect  the 
official  views  or  policy  of  the  Department  of  Transportation.  This 
report  does  not  constitute  a standard,  specification  or  regulation. 


pst*5scs  ?*?  s»? 


wqF353JJ?w* ,,  ,«.5^-M',“'--  '•• 


. . „ »^vi, A»vrj^^t^’‘,,n^WrKe  »■“* 


PREFACE 


The  advice  and  assistance  of  associates  at  Booz,  Alien  Applied 
Research  was  sincerely  appreciated:  Mr.  John  Wing.  Mr.  David  Weiss. 
Dr.  Russ  Kevala,  and  Mrs.  Sharon  Burns  Holcombe. 


The  authors  would  like  to  thank  all  those  individuals  and  organiza- 
tions who  provided  data  for  this  study  and  for  our  predecessor  study  on 
glass  upon  which  this  one  is  based: 

National  Bureau  of  Standards: 

Dr.  Edward  B.  Magrab 

Mr.  E.  V.  Lyndecker 

Mr.  James  Gross 

Mr.  James  Cliston 

Mr.  Webster  Capps 

Dr..  Sheldon  M.  Wiederhorn 

Mr.  James  R.  Clifton 

National  Concrete  Masonry  Association: 

Mr.  Thomas  Redmond 

National  Gypsum  Company: 

Mr.  David  Bicri 

Gypsum  Association: 

Mr.  Herbert  B.  Carlson 
Mr.  C.  M.  Hedlund 

U.  S.  Gypsum  Company 

Mr.  J.  M.  Sommerfield 

Brick  Institute  of  America: 

Mr.  Alan  Yorkdale 

National  Lime  Association: 

Mr.  Kenneth  Gutschick 

NASA  Langley: 

Dr.  Harvey  H.  Hubbard 
Mr.  Domenic  J.  Maglieri 
Mr.  Harry  W.  Carlson 
Mr.  Arnold  W.  Mueller 


. ***» 


John  A.  Blume  and  Associates: 
Mr.  Lloyd  Lee 

Wyle  Laboratories: 

Mr.  Robert  White 


Libbey-Owens-Ford  Company 

PPG  Industries: 

Mr.  Robert  W.  McKinley 

Corning  Glass  Works: 

Dr.  John  S.  McCartney 
Dr.  John  Munier 

Alfred  University: 

Dr.  Charles  H.  Greene 
Dr.  Harry  Stevens 

Air  Force  Office  of  Research: 
Mr.  Mathew  J.  Kerper 

Naval  Research  Laboratory: 
Dr,  Steve  Freiman 

Virginia  Polytechnic  Institute: 
Dr.  Robert  Heller 

University  of  Southhampton: 

Dr.  G.  M.  Lilley 

Royal  Aircraft  Establishment: 
Mr.  David  R.  B.  Webb 


iv 


s 

i 


r 


i 


TABLE  OF  CONTENTS 


INTRODUCTION 


SONIC  BOOM  I DA  DING 


Page 

Number 


1. 

Nominal  Overpressure,  pQ 

3 

2. 

Free  Field  Overpressure,  pf 

6 

3. 

External  Overpressure.  pe 

6 

4. 

Internal  Overpressure,  pj 

7 

5. 

Static  Stress 

9 

6. 

Dynamic  Amplification  Factor  (DAF) 

9 

7. 

Sonic  Boom  Stress  Analysis 

11 

III. 

SONIC  BOOM  EXPERIMENTS 

16 

1. 

Cklahoma  City  Tests 

16 

2. 

White  Sands  Tests 

28 

3. 

Edwards  Air  Force  Base  Tests 

36 

IV. 

RESPONSE  PROBABILITY  DENSITY 
FUNCTION  TECHNIQUE 

45 

1. 

Probability  Density  Functions 

45 

2. 

Sonic  Boom  Stress 

54 

3. 

Material  Strength 

57 

4. 

Probability  of  Breakage 

59 

V. 

GLASS 

63 

VI. 

PLASTER 

86 

v 


w w m u ^ ^ «?  * * twwwttp  7 -'^n  r>*-[  't'^J^,,T 


12. 


13. 


14. 


15. 

16. 

17. 

18. 

19. 

20. 


21. 


22. 


23. 

24. 


Probability  of  Equaling  or  Exceeding  a Given 
Value  of  the  Ratio  of  Measured  to  Calculated 
Overpressures  for  Airplane  B at  Station  3 


Probability  of  Equaling  or  Exceeding  a Given 
Value  of  the  Ratio  of  Measured  to  Calculated 
Overpressures  for  Airplane  B at  Station  4 


Instrument  Locations  for  White  Sands  Test 
Test  Structure  W4 


One-Story  Test  Structure  (E-l)  EAFB 
Two-Story  Test  Structure  (E-l)  EAFB 
Example  of  a Probability  Density  Function 
Probability  Density  Function  of  Stress 
Probability  Density  Function  of  Breaking  Strength 


Probability  Density  Function  of  the  Effective 
Factor  of  Safety  for  the  Lognormal  Case 


Experimental  Probability  Density  Functions  of 
Logarithm  of  External  Fressure/Nominal 
Pressure.  Station  3 


Development  of  i icoability  Density  Function 
of  log10 


Stress-Time  Characteristics  of  Glass 


Comparison  of  Breaking  Pressure  Distributions 
of  Used  and  New  Glass 


25.  Model  of  the  Probability  Density  Function  of 
the  Logarithm  of  the  Breaking  Strength  of 
the  Glass  Population 


vm 


Page 

Number 


26 


27 


30 

37 

38 
46 

49 

50 


53 


58 


61 

65 


67 


69 


; £ 


i i 

-e 

\ i. 


1 5 

!i 


5 1 


! \ 


% a 

\ j 


Page 

Number 


Breakage  Probabilities  for  Glass  in  Good  Condition 
for  All  Flight  Paths  Equally  Likely 

Breakage  Probabilities  for  Glass  in  Good  Condition 
for  Head-Cn  Flight  Paths 

Variation  in  Mechanical  Properties  of  Gypsum 
Plasters  With  Consistency 

Stress  Strain  Plot  for  Plaster  (80  Parts  Water) 


Breakage  Probabilities  for  Plaster  for  All  Flight 
Paths  Equally  Likely 

Breakage  Probabilities  for  Brie -a -Brae  for 
All  Flight  Paths  Equally  Likely 

Breakage  Probabilities  for  Free-Standing  4-Inch 
Brick  W'alls  for  All  Flight  Paths  Equally  Likely 

Ranges  of  Breakage  Probabilities  for  Glass,  Plaster. 
Bric-a-Brac.  and  Free-Standing  Brick  Wails  of 
High-Bond  Mortar 


INDEX  OF  TABLES 

Typical  Oklahoma  City  Data 

White  Sands  Strain  Readings,  Rosette  R-l, 
November  27,  1964 


White  Sands  Strain  Readings,  Rosette  R-2,  j ■' 

November  27,  1964  32  1 1 


White  Sands  Strain  Readings,  Rosette  R-l, 

November  28,  1964  33 


White  Sands  Strain  Readings,  Rosette  R-2, 
November  28  1964 

Angie  of  Incidence  on  Structure  W4 

External  Pressures  on  Structure  W4  at  White  Sands 

Summary  of  Instrumented  Windows 

Summary  of  Edwards  Air  Force  Base  Data 

Mean  of  the  Logarithm  of  the  Effective  Factor  of 
Safety  for  All  Angles  of  Incidence  Combined 

Variance  of  the  Logarithm  of  the  Effective  Factor 
of  Safety  for  All  Angles  of  Incidence  Comb.-*ed 

Mean  of  the  Logarithm  of  the  Effective  Factor 
of  Safety  for  0°  Angle  of  Incidence 

Variance  of  the  Logarithm  of  the  Effective  Factor 
of  Safety  for  0°  Angle  of  Incidence 


Page 

Number 


Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  Ail  Flight  Paths 

Equally  Likely  at  a Nominal  Overpressure 

of i psf  77 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  Head-On  Flight  Path 

at  a Nominal  Overpressure  of  1 psf  78 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  All  Flight  Paths 

Equally  Likely  at  a Nominal  Overpressure 

of  2 psf  79 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  Head-On  Flight  Path 

at  a Nominal  Overpressure  of  2 psf  80 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  Head-On  Flight  Path 

at  a Nominal  Overpressure  of  4 psf  81 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  fo t Head-On  Flight  Path 

at  a Nominal  Overpressure  of  20  psf  82 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  Head-On  Flight  Path 

at  a Nominal  Overpressure  of  40  psf  83 

Probability  of  Glass  Breakage  Calculated  for 

Model  Window  Population  for  Head-On  Flight  Path 

at  a Nominal  Overpressure  of  100  psf  84 

Strength  of  Basecoat  Plaster  for  Mixtures 

With  Various  Aggregates  90 

Plaster  Breakage  at  White  Sands  93 

Static  Failure  Pressures  of  Plaster  Configurations  95 


Page 

Number 


25. 

Damage  Claim  Tabulation  for  Edwards  Air  Force 
Base  Sonic  Boom  Tests 

103 

26. 

Approximate  Ranges  in  Compressive  Strength  for 
Various  Mortar  Proportions  and  Types 

107 

27. 

Conditions  of  Failure  of  Peak  Overpressure- 
Sensitive  Elements 

109 

28. 

Static  Failure  Pressures  of  Brick 
Wall  Configurations 

110 

29. 

Measures  to  Reduce  Sonic  Boom  Structural 
Breakage  Rate 

117 

-JTTTT-?.' 


I.  INTRODUCTION 

There  is  a considerable  body  of  data  on  sonic  boom  structural 
damage  effects  from  various  experiments  that  have  been  conducted 
over  the  past  fifteen  years.  Experiments  have  shown  that  there  is  a 
small  but  finite  probability  of  sonic  booms  breaking  especially  sus- 
ceptible items  such  as  windows,  plaster,  and  bric-a-brac.  The 
probability  of  breakage  was  consistently  found  to  increase  with  higher 
overpressure  sonic  booms.  The  results  of  the  experiments  have  in 
general  agreed  with  the  Air  Force's  experience  in  reimbursing 
damage  claims  from  the  public  from  supersonic  overflights.  The 
claims  most  frequently  substantiated  and  reimbursed  have  been  for 
damage  to  glass,  plaster,  and  bric-a-brac. 

Window  glass  is  generally  the  most  susceptible  element  in  a 
house  to  sonic  boom  damage.  Our  previous  work  [1,2, 3, 4]  has 
determined  the  probability  of  breakage  of  various  window  configura- 
tions as  a function  of  sonic  boom  overpressure.  It  was  calculated 
from  our  statistical  model  that  for  a typical  window  population  there 
would  be  1.1  breaks  per  million  panes  in  good  condition  boomed  at 
a nominal  overpressure  of  1 pound  per  square  foot  (psf).  This 
est;mate  agreed  well  with  sonic  boom  claims  data. 
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The  statistical  modeling  method  which  was  used  in  the  previous 
studies  of  glass  breakage  is  a method  we  call  "the  response  probability 
density  function  technique.  " This  technique,  which  will  be  described 
in  detail  in  another  chapter,  combines  the  probability  density  function 
(pdf)  of  the  sonic  boom  stress  with  the  pdf  of  the  material  strength. 

The  resultant  pdf  readily  yields  the  probability  of  structural  failure. 

Since  the  response  pdf  technique  was  first  applied  to  sonic 
booms  on  glass  in  1972-73  we  have  made  substantial  progress  in 
extending  it  to  other  structural  response  applications.  It  was  applied 
to  estimating  the  probability  of  air  terminal  window  breakage  from 
the  noise  from  taxiing  aircraft  in  1974  [5] . Then  in  1975  it  was 
applied  in  the  Concorde  environmental  impact  statement  to  determine 
the  breakage  probabilities  of  windows,  plaster  ceilings,  stone  bridges, 
and  brick  chimneys  to  noise  from  subsonic  overflights.  On  the  basis 
of  the  success  of  the  technique  for  these  cases,  this  study  extends  the 
method  to  piv  ''^cr.  c?  sonic  boom  breakage  for  plaster,  bric-a-brac, 
and  brick.  Thus,  with  the  glass  previously  accounted  for  and  the 
three  items  considered  in  this  study,  there  is  now  a comprehensive 
treatment  for  predicting  the  probability  of  breakage  for  the  vast 
majority  of  sonic  boom  structural  damage  cases. 
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II.  SONIC  BOOM  LOADING 

The  basic  problem  in  sonic  boom  breakage  prediction  is  one 
o * finding  the  stress  created  by  a transient  overpressure.  The  over- 
pressure caused  by  a sonic  boom  is  usually  an  N-wave  (see  Fig- 
ure 1).  Those  small  variations  in  the  wave  shape  that  make  some 
individual  waves  more  peaked  or  more  rounded  than  the  ideal  N-wave 
have  little  influence  on  the  results  derived  here.  These  variations 
are  accounted  for  in  the  fact  that  the  present  statistical  model's 
parameters  were  derived  from  actual  field  data  which  represented 
the  typical  waveforms  encountered.  The  paragraphs  which  follow 
review  the  definitions  of  the  parameters  shown  in  Figure  2. 

1.  NOMINAL  OVERPRESSURE,  pQ 

Adequate  methods  presently  exist  for  calculating  a predicted 
peak  sonic  boom  overpressure  value  based  on  aircraft's  altitude, 
speed,  weight,  and  dimensions  (7J.  This  calculated  value  is  called 
the  nominal  overpressure,  po.  Its  value  is  generally  in  the  range 
of  1 psf  to  4 psf  for  typical  supersonic  aircraft  operations.  It  is 
desired  to  predict  P,  the  probability  of  a window  being  broken  by  a 
given  boom,  as  a function  of  pQ. 
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Figure  2.  Evolution  of  Sonic  Boom  Load 


2.  FREE  FIELD  OVERPRESSURE,  pf 

The  sonic  boom  suffers  various  attenuation  and  focusing  effects 
in  passing  through  the  atmosphere  to  the  ground.  The  random  inhomo- 
geneities in  the  atmosphere  cause  a spread  in  the  range  of  values  of 
the  peak  free  field  overpressure  (p^)  that  is  measured  on  the  ground. 

By  convention,  the  value  of  pf  is  measured  with  a microphone  mounted 
on  the  ground,  usually  on  a plywood  board.  This  method  of  measure- 
ment results  in  a pressure  doubling  effect  from  the  ground  reflection 
being  included  in  the  value  of  pf.  Thus  the  measured  value,  pf,  is 
not  really  a free  field  pressure  as  the  term  is  commonly  used  in 
acoustics,  but  it  is  conventionally  called  this  in  the  sonic  boom  com- 
munity. Similarly,  the  calculation  of  pD,  the  nominal  pressure,  in- 
cludes a conventional  value  of  the  ground  reflection  coefficient  of 
1.9.  In  general,  the  mean  value  of  pf  will  be  very  close  to  p0,  but 
individual  values  of  pf  will  exhibit  considerable  spread.  This  random 
effect  results  from  the  atmospheric  inhomogeneities  and  from  dif- 
ferences in  the  aircraft  operating  parameters  from  those  used  for 
calculating  po- 

3.  EXTERNAL  OVERPRESSURE,  pe 

When  a sonic  boom  having  a free  field  peak  overpressure  pf 
impinges  on  a structure  it  creates  an  external  overpressure  waveform 
having  a peak  value  pe.  The  value  of  pe  will  depend  on  the  angle  the 
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incoming  wave  makes  with  the  structure,  the  reflection  coefficient 
of  the  structure  surface,  and  whether  the  wave  combines  with  other 
waves  reflected  off  other  structures  or  the  ground.  For  instance,  it 
has  generally  been  observed  In  experiments  f8,  9)  that  the  side  of 
a house  exposed  to  a head-on  sonic  boom  wave  exhibits  a much 
higher  overpressure  than  the  opposite  side.  There  have  been  some 
attempts  [10,  11, 12]  to  obtain  analytical  solutions  for  pc  in  terms  of 
p^  by  modeling  the  acoustical  scattering.  These  attempts,  however, 
involved  complicated  computer  solutions  which  were  applicable  only 
to  a single  structure  geometry  and  wave  angle.  For  practical  engi- 
neering situations,  the  relation  between  pf  and  pe  would  seem  to  be 
best  treated  by  statistical  methods. 

4.  INTERNAL  OVERPRESSURE,  pj 

TIk  waveform  of  the  overpressure  inside  a structure  depends 
on  the  transmissibility  of  the  structure  to  sound  as  a function  of  its 
geometry.  The  solution  for  the  internal  waveform  is  quite  compli- 
cated and  differs  for  each  structure  [13, 14].  It  is  generally  observed 
[13, 14]  that  the  internal  overpressure  waveform  is  sinusoidal  with 
a period  the  same  as  that  of  the  N-wave  and  a peak  value  pj  which  is 
some  fraction  of  the  external  peak  overpressure  pe.  Subtracting  the 
waveform  of  the  internal  pressure  from  that  of  the  external  pressure 
gives  the  net  loading  waveform  acting  on  the  window  (see  the  top  of 
Figure  3). 
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STATIC  STRESS 
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The  analysis  of  the  loading  of  a structure  by  a sonic  boom  has 
been  a problem  that  has  attracted  much  interest  in  the  past  fifteen 
years  f9,  14-21],  Generally  these  analyses  have  included  a compari- 
son with  loading  under  uniform  static  pressure.  This  is  because  most 
structural  elements  are  designed  to  static  wind  loading,  and  hence 
static  pressure  design  data  are  much  more  readily  available  than  some 
boom  loading  data. 

A window  configuration,  for  instance,  is  usually  assumed  to  be 
i rectangular  plate  simply  supported  at  all  four  edges;  which  has  been 
found  to  agree  with  experimental  studies  of  static  window  loading  f 22] . 
For  small  deflections,  the  stress  is  a linear  function  of  the  static 
pressure  and  is  given  by  the  well-known  formula  [23] 


2h2(a2+b2) 

where  ag  is  the  extreme  fiber  stress  at  the  center  of  the  plate,  ps  is 
the  static  pressure,  a is  the  longer  side  of  the  rectangle,  b is  the 
shorter  side,  and  h is  the  thickness  of  the  plate. 


6,  DYNAMIC  AMPLIFICATION  FACTOR  (DAF) 

Since  windows  are  usually  designed  to  static  load  requirements, 
it  is  desirable  to  compare  sonic  boom  loads  to  such  static  loads.  This 
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is  accomplished  by  finding  the  dynamic  amplification  factor  (DAF) 
defined  by 

DAF=~L.  (2) 

a d 

The  numerator  of  the  DAF  is  the  maximum  stress  in  the  plate  and 

the  denominator  is  the  dynamic  stress  defined  by 

2 2 
p a b 

ad  ■>  - Pe  F (3) 

2h  (a  +b“) 

where  F is  a stress  factor.  Note  that  a,  is  tne  stress  that  would  be 

d 

produced  by  a static  pressure  equal  in  magnitude  to  the  pea',  external 
pressure  p^. 

Throughout  the  literature  there  are  various  definitions  of  DAF. 
Some  definitions  use  a ratio  of  deflections  rather  than  stresses,  and 
others  use  the  peak  value  of  the  net  loading  waveform  pnet  rather 
than  pe,  the  peak  value  of  the  external  pressure  waveform.  Since 
the  problem  at  hand  is  one  of  rupture  and  not  deflection  the  current 
definition  was  deemed  more  appropriate.  In  the  definition  of 
of  (3)  it  was  decided  to  use  p0  rather  than  pnet  because  (a)  pe  is 
more  easily  measured,  (b)  pe  is  very  nearly  equal  to  pnet  as  can  be 
seen  from  the  top  of  Figure  3,  (c)  the  available  overflight  data  base 
contains  only  values  of  pe#  and  (d)  it  was  found  by  statistical  analysis 
that  the  internal  pressure  does  not  correlate  well  with  window  stress. 


7. 


SONIC  BOOM  STRESS  ANALYSIS 


According  to  theory,  when  a thin  plate  such  as  a window  is  loaded 
to  the  extent  that  its  deflection  exceeds  one  half  its  thickness,  nonlinear 
effects  start  to  come  into  play.  This  results  because  membrane  stresses 
become  comparable  in  magnitude  to  the  bending  stresses.  Because  of 
this  effect,  the  stress  will  be  somewhat  lower  than  that  predicted  by 
Equation  1 once  the  deflection  starts  to  exceed  one  half  the  thickness. 

The  nonlinear  differential  equations  for  large  deflections  have  been 
solved  in  recent  years  [24,25]  and  the  resultj  have  beer,  presented  as  a 
series  of  graphs  for  rectangles  having  various  aspect  ratios  [26]. 


Unfortunately  the  stresses  in  glass  plates  correlate  very  poorly 
with  large  deflection  theory  [14].  Although  the  center  deflections  of  the 
plates  follow  the  nonlinear  theory  fairly  well  the  stresses  do  not.  Thus 
what  some  authors  [14-16]  have  done  is  to  calculate  DAF's  on  the  basis 
of  deflection  using  dynamic  analysis  with  the  nonlinear  theory.  It  can 
be  seen  from  Figures  3 and  4 from  the  linear  theory  and  Figure  5 from 
the  nonlinear  theory  that  the  envelopes  of  the  DAF  curves  seem  to  range 
between  1 and  2 whether  the  linear  or  nonlinear  theory  is  used.  The 
curves  in  Figures  3 and  4 were  developed  by  Seaman  [14]  using  the 
conventional  approach  of  modeling  the  window  as  a linear  single  degree 
of  freedom  system,  finding  its  impulse  response,  and  convolving  with 
the  excitation  to  find  the  maximum  excursion.  The  slight  difference 
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between  the  curves  in  Figures  3 and  4 results  from  the  fact  that  Figure  3 

Pi 

uses  an  internal  pressure  ratio  q = — = . 25  and  Figure  4 has  q - . 5. 

© 

Since  these  figures  resulted  from  linear  analysis  the  deflection  DAF's 
shown  are  numerically  equal  to  stress  DAF's.  Seaman  then  went 
through  a large  deflection  analysis  which  resulted  in  the  curve  in 
Figure  5.  The  other  references  cited  [15-21]  all  displayed  results 
similar  to  Seaman's.  These  three  curves  are  introduced  here  to  point 
out  two  striking  similarities  which  hold  regardless  of  whether  deflection 
is  assumed  iinecr  or  nonlinear  and  no  matter  how  great  the  value  of 


Most  of  the  values  of  the  DAF  are  between  1 and  2. 

The  values  of  the  DAF  below  1 occur  only  for  small  values 
of  the  period  rati  for  cases  where  a structural  element 
with  a low  resonant  frequency  fQ  is  excited  by  a short 
duration  boom  with  a small  value  of  T. 


The  second  point  above  emphasizes  that  the  structural  element's 
frequency  can  be  important  in  some  cases.  Assuming  linear  deflection 
and  simple  support  around  all  four  edges  of  a plate,  its  natural  fre- 


quency is  [27) 


f -ill  / y ... 

" 2 v 12p(l-y2) 


ih  ■&) 


where  Y is  Young's  modulus,  p is  the  mass  density,  and  v is  Poisson's 
ratio.  A case  of  a low  period  ratio  can  occur  when  an  F-104  aircraft 
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with  T = . 067  sec  booms  a large  storefront  window  with  f0  » 4Hz.  In 
this  case,  the  window  suffers  minimal  deflection  because  the  excita- 
tion is  considerably  above  the  window's  natural  frequency. 

In  the  present  study  the  linear  stress  theory  of  Equation  (1)  is 
used  exclusively.  The  reasons  for  disregarding  the  large  deflection 
theory  here  are  the  following: 

1.  The  DAF  measurements  used  in  the  present  study  were 
at  small  deflections  in  the  linear  range. 

2.  The  large  deflection  theory  gives  poor  results  for  pre- 
dicting breaking  stress. 

3.  The  coefficient  relating  stress  to  pressure  in  Equation  (1) 
will  be  seen  to  cancel  out  later  in  the  analysis  anyway. 

In  effect  we  will  be  working  only  with  boom  pressure  and 
breaking  pressure,  assuming  the  relation  between  stress 
and  pressure  to  be  linear  between  the  sonic  boom  test 
pressures  and  the  structural  element  breaking  pressures. 
This  allows  the  use  of  linear  transformations. 

Having  surveyed  the  problem  of  sonic  boom  loading,  the  dis- 
cussion will  now  review  the  data  from  sonic  boom  structural  response 
experiments. 
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III.  SONIC  BOOM  EXPERIMENTS 


A large  data  base  consisting  of  over  70,  000  readings  from 
various  experiments  of  the  past  twelve  years  has  been  utilized  in 
this  study.  The  experiments  which  generated  this  data  were  the 
following: 

1.  Oklahoma  City  Sonic  Boom  Tests  performed  by  FAA  and 
NASA  in  1964  [28]  . 


2.  White  Sands  Sonic  Boom  Structural  Reaction  Tests  per- 
formed by  FAA  in  1965  [29,30] . 

3.  Edwards  Air  Force  Base  Sonic  Boom  Tests  performed 
by  FAA  and  NASA  in  1966  [31,32]  . 


Each  of  these  experiments  will  now  be  described,  together  with  the 
contribution  of  its  data  to  the  present  study. 


1.  OKLAHOMA  CITY  TESTS 

The  sonic  boom  tests  conducted  in  the  Oklahoma  City  area  from 
April  to  July  1964  were  primarily  to  determine  community  response. 
In  addition  to  surveying  how  acceptable  the  booms  were  to  the  public, 
the  tests  also  included  measurements  of  ground  overpressure.  The 
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overpressures  measured  were  values  of  p^  , the  peak  free  field 
overpressure.  Over  3500  sonic  boom  signatures  were  recorded  during 
this  test  series.  These  signatures  form  the  basis  of  the  present  analysis 
of  the  relation  between  peak  nominal  pressure  pQ  and  peak  free  field 
pressure  p^. . 


The  aircraft  used  to  produce  the  Oklahoma  City  sonic  booms  were 
the  F-104,  F-101,  F-106,  and  B-58  shown  in  Figure  6.  These  were 
planes  whose  characteristics  had  been  thoroughly  analyzed  by  NASA, 
and  they  were  able  to  calculate  a nominal  overpressure  p^  for  each 
flight.  The  values  of  pQ  ran  from  .64  psf  to  2.  17  psf  for  the  F-101 
and  F-104  flights.  Only  the  F-101  and  F-104  data  are  analyzed  here  since 
there  were  an  insufficient  number  of  readings  from  F-1061  s and  B-58's 
to  be  statistically  meaningful. 


The  flight  track  for  the  test  series  is  shown  in  Figure  7.  Note 
that  there  were  three  ground  measurement  stations  at  various  distances 
from  the  flight  track:  Station  1,  Station  3,  and  Station  4.  A test  building 
was  located  at  each  station  with  a microphone  shock  mounted  inside. 
These  inside  microphones  which  measured  pi  were  at  the  center 
of  the  room  5 ft  from  the  floor  level.  Outside  microphones,  which 
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measured  pf  were  shock-mounted  at  ground  level  on  the  surface  of  a 

plywood  reflecting  board.  NASA  obtained  strip  charts  of  internal  and 

free  field  overpressure  for  each  boom  at  each  station.  From  these 

strip  charts  NASA  determined  the  peak  values  p,  and  p and  the 

f 1 

positive  duration  T/2  and  positive  impulse  1/2  for  each  boom.  A 
representative  page  of  Oklahoma  City  data  is  shown  in  Table  1.  Note 
that  ApQ , Apj , At^pos,  and  I^pos  of  the  column  headings  are  p^.  , 
p.  , T/2  , and  1/2  , respectively,  in  the  present  notation.  The 
nominal  overpressure  values  for  each  station  have  been  added  at  the 
bottom  of  each  station  set  for  comparison. 

In  its  analysis  of  the  data  NASA  calculated  the  ratio  p /p  for 

f o 

each  boom  measured.  Histograms  of  this  ratio  were  then  plotted  by 
NASA  for  each  station's  readings  for  the  F-101  and  F-104.  These  six 
histograms  are  shown  in  Figures  8-13  In  the  present  notation 

APomeas/ APocalc  of  the  abscissas  is  Pf/PQ.  Each  figure  also 
includes  a cumulative  probability  plot  of  the  data  on  a lognormal  scale. 

Note  from  the  straight  lines  shown  for  comparison  that  the  data  do 
appear  to  be  distributed  lognormal.  As  will  be  discussed  later 
this  lognormal  probability  density  leads  to  a simplification  of  the 
analysis. 
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Table  1.  Typical  Oklahoma  City  Data 
Summary  of  Sonic -Boom  Data  for  Various  Flights  of  Airplane  A 
for  an  Altitude  of  32,  000  Feet 
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•Fot  thete  two  leading!  nominal  ovetpiesjum  wtie:  Station  l,p0  » 1.27;  Station  3,  p0  » 1.09;  Station  4,  p0  *■  .81. 
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Figure  8.  Probability  of  Equaling  or  Exceeding  a Given  Value  of  the 
Ratio  of  Measured  to  Calculated  Overpressures  for 
Airplane  A at  Station  1 
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Figure  10.  Probability  of  Equaling  or  Exceeding  a Given  Value  of  the 
Ratio  of  Measured  to  Calculated  Overpressures  for 
Airplane  A at  Station  4 
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•’igure  11.  Probability  of  Equaling  or  Exceeding  a Given  Value  of  the 
Ratio  of  Measured  to  Calculated  Overpressures  for 
Airplane  B ?*.  Station  1 
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Figure  12.  Probability  of  Equaling  or  Exceeding  a Given  Value  of  the 
Ratio  of  Measured  to  Calculated  Overpressures  for 
Airplane  B at  Station  3 
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conducted  by  FAA  at  the  White  Sands  test  range.  Se . eral  structures 
were  equipped  with  microphones  on  various  external  walls  to  measure 


external  pressure  pg  . In  addition  there  were  also  microphones  outside 
on  the  ground  for  measuring  and  inside  at  the  centers  of  rooms 
for  measuring  Pj  . For  measuring  the  structural  reactions  there 
were  strain  gauges,  scratch  gauges,  and  accelerometers  attached  to 
windows,  walls,  and  roofs. 

The  aircraft  used  in  the  White  Sands  tests  were  F-104^  and  B-58’s 
as  shown  in  Figures  6a  and  6d.  These  aircraft  were  flown  over  the 
Uninhabited  test  structures  at  much  lower  altitude  than  for  normal 
operation*  to  deliberately  create  high  overpressures.  While  typical 
sonic  booms  from  nc  -mal  aircraft  operations  are  of  the  order  of  1-4  psf, 
those  in  this  test  ceries  were  often  in  the  10-20  psf  range.  These  high 
o ■ ,*pre8surea  were  necessary  to  cause  the  large  structural  responses 
necessary  to  explore  damage  conditions. 
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The  particular  structure  which  provided  the  data  used  in  the  present 
study  was  designated  W4.  The  instrumentation  on  this  structure,  shown 
in  Figure  14,  consisted  of  strain  gauges  (labeled  R),  microphones 
(labeled  M),  velocimeters  (labeled  V),  and  accelerometers  (labeled  A). 

In  addition  to  the  instrumentation  shown  there  were  also  an  interior 
microphone  M-7  and  an  exterior  ground-mounted  free  field  microphone 
M-16.  The  readings  from  the  instrumentation  which  were  utilized  in 
the  present  study  are  shown  in  Tables  2-5  [33]  . In  each  case  the  table 
values  of  strain  and  pressure  are  the  peak  values  for  the  tape  recorded 
sequence.  Note  that  the  column  headed  "vector"  has  the  letters  A-H. 
These  correspond  to  eight  possible  angles  of  sound  wave  incidence  in 
45°  increments  as  keyed  to  Table  6. 


In  addition  to  the  foregoing  strain  gauge  data  from  building  W -4, 

there  was  also  another  series  of  484  external  pressure  reading  sets 

from  this  same  building.  These  readings  formed  the  basis  of  the  present 

statistical  analysis  of  the  effect  of  the  angle  of  incidence.  The  first  30 

of  these  sets  of  readings  are  shown  in  Table  7 as  printed  out  by  the 

computer,  including  some  transformed  values.  The  fourth  column 

contains  the  cosines  of  the  angles  between  the  incoming  sound  rays  and 

the  normal  to  the  surface.  Note  from  the  values  of  the  cosines  that  as 

in  Table  6 all  angles  have  been  approximated  as  the  nearest  45°  multiple. 

Columns  5 and  7 contain  the  ratio  p /p.  for  the  north  and  south  walls, 
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Microphone  M-7  located  in  middle  of  room 
interior;  3 feet  behind  R-l  and  R-2. 


Figure  14.  Instrument  Locations  for  White  Sands  Test  Structure  W4 
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Bun 

1- 145 

2- 146 

3- 147 

4- 148 

5- 149 

6- 150 

7- 151 

8- 152 

9- 153 

10- 154 

11- 155 

12- 156 

13- 157 

14- 158 

15- 159 

16- 160 

17- 161 

18- 162 

19- 163 

20- 164 

21- 165 

22- 166 

23- 167 

24- 168 

25- 169 

26- 170 

27- 171 

28- 172 

29- 173 

30- 174 


Table  2.  White  Sands  Strain  Readings  Rosette  R-l, 
November  27,  1964 


Mach 

Alt. 

Streln  u 

In. /In. 

Recorded  Pressures  - osf 

x 1.000 

Vector 

a 

b 

c 

M-l 

M-7 

M-16 

1.40 

13.5 

D 

57 

35 

NR 

NR 

NR 

NR 

1.43 

14.0 

H 

54 

29 

NR 

3.1 

1.9 

3.5 

1.32 

13.5 

D 

NR 

NR 

NR 

2.4 

2.7 

3.9 

1.30 

14.0 

E 

40 

26 

NR 

1.4 

0.9 

3.9 

1.30 

14.0 

c 

80 

48 

NR 

NR 

NR 

NR 

1.30 

14.0 

E 

69 

42 

NR 

1.9 

0.9 

3.9 

1.39 

15.0 

F 

40 

26 

NR 

1.8 

0.9 

3.4 

1.34 

13.0 

A 

66 

42 

NR 

5.5 

1.1 

3.4 

1.38 

14.0 

F 

46 

29 

NR 

1.9 

1.3 

3.1 

1.40 

14.0 

Q 

F2 

29 

NR 

2.2 

0.9 

3.9 

1.40 

14.3 

C 

43 

29 

NR 

1.8 

1.8 

3.9 

1.38 

14.0 

G 

NR 

NR 

NR 

NR 

NR 

NR 

1.30 

14.0 

K 

86 

48 

NR 

3.6 

1.4 

2.0 

1.36 

14.0 

D 

66 

35 

NR 

2.5 

2.2 

3.4 

1.36 

14.0 

H 

57 

42 

NR 

3.2 

1.4 

1.9 

1.25 

13.5 

A 

66 

35 

NR 

5.0 

2.2 

7.4 

1.30 

14.0 

F 

63 

36 

NR 

1.7 

2.2 

3.7 

1.30 

14.0 

A 

89 

52 

NR 

3.2 

1.2 

2.7 

1.20 

14.0 

B 

100 

58 

NR 

3.2 

1.9 

2.5 

1.38 

is.n 

E 

60 

36 

NR 

1.5 

1.0 

1.9 

1.34 

14.0 

B 

83 

45 

NR 

4.3 

2.0 

5.0 

1.45 

16.0 

C 

77 

42 

NR 

NR 

NR 

NR 

1.39 

14.3 

6 

72 

39 

NR. 

NR 

NR 

NR 

1.45 

13.9 

C 

80 

45 

NR 

NR 

NR 

NR 

1.30 

14.0 

0 

80 

45 

NR 

4.1 

2.7 

3.9 

1.30 

14.0 

H 

117 

65 

NR 

3.2 

1.4 

5.3 

1.40 

14.0 

0 

69 

36 

NR 

2.3 

1.8 

3.9 

1.35 

14.0 

E • 

34 

19 

NR 

1.9 

0.8 

2.0 

1.30 

14.0 

B 

72 

42 

NR 

4.5 

1.4 

3.2 

1.33 

14.0 

E 

23 

13 

NR 

1.7 

0.8 

4.3 

31 
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Table  3.  White  Sands  Strain  Readings  Rosette  R-2, 
November  27,  1964 


Streln  u In. /In.  Recorded  Pressures  - psf 


Run 

Hech 

HI  »• 

x 1.000 

Vector 

« 

b 

c 

M-2 

M-7 

M-16 

1-145 

1.40 

13.5 

D 

19 

35 

28 

NR 

NR 

NR 

2-146 

1.43 

14.0 

H 

11 

9 

NR 

4.6 

1.9 

3.5 

4-145 

1.30 

14.0 

E 

7 

9 

6 

1.9 

0.9 

3.9 

5-149 

1.30 

14.0 

B 

30 

39 

25 

NR 

NR 

NR 

6-150 

1.30 

14.0 

E 

11 

13 

6 

2.3 

0.9 

3.9 

7-151 

1.39 

15.0 

F 

11 

13 

9 

3.5 

0.9 

3.4 

8-152 

1.34 

13.0 

A 

19 

22 

31 

4.5 

1.1 

3.4 

9-153 

1.38 

14.0 

F 

11 

26 

19 

4.6 

1.3 

3.1 

10-154 

1.40 

14.0 

G 

22 

30 

28 

1.5 

0.9 

3.9 

11-155 

1.40 

14.3 

C 

7 

13 

9 

1.9 

1.8 

3.9 

13-157 

'..30 

14.0 

H 

33 

44 

46 

3.3 

1.4 

2.0 

’4-158 

1.36 

14.0 

D 

22 

22 

22 

3.6 

2.2 

3.4 

15-159 

1.36 

14.0 

. H 

19 

35 

31 

2.5 

1.4 

.9 

16-160 

1.25 

13.5 

A 

30 

39 

38 

6.5 

2.2 

7.4 

17-161 

1.30 

14.0 

F 

37 

44 

34 

2.6 

2.2 

3.7 

18-162 

1.30 

14.0 

A 

19 

30 

25 

4.5 

1.2 

2.7 

19-163 

1.20 

14.0 

B 

30 

44 

41 

2.4 

1.9 

2.5 

20-164 

1.38 

15.0 

E 

19 

22 

25 

3.2 

1.0 

1.9 

21-165 

1.34 

14.0 

B 

22 

35 

28 

2.3 

2.0 

5.0 

22-166 

1.45 

16.0 

C 

11 

13 

13 

NR 

NR 

NR 

23-167 

1.39 

14.3 

G 

19 

26 

28 

NR 

NR 

NR 

24-168 

1.45 

13.9 

C 

7 

9 

13 

NR 

NR 

NR 

25-169 

1.30 

14.0 

0 

11 

13 

9 

2.9 

2.7 

3.9 

26-170 

1.30 

14.0 

H 

26 

49 

44 

5.6 

1.4 

5.3 

27-171 

1.40 

14.0 

D 

7 

9 

13 

2.0 

1.8 

3.9 

28-172 

1.35 

14.0 

E 

7 

9 

9 

1.8 

0.8 

2.0 

29-173 

1.30 

14.0 

B 

30 

39 

34 

6.0 

1.4 

3.2 

30-174 

1.33 

14.0 

E 

7 

9 

9 

2.3 

0.8 

4.3 
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Table  4.  White  Sands  Strain  Readings  Rosette  R-l,  i- 

November  28,  1964  ii 

li 


Rub 

Mach 

Alt. 
x 1.000 

Victor 

Strain  u 

a b 

In. /In. 
c 

Recorded  Pressures  - psf 
M-l  M-7  N-16 

) 

2-176 

1.30 

11.9 

A 

109 

61 

NR 

9.1 

3.6 

5.9 

i 

3-177 

1.25 

11.7 

F 

52 

29 

NR 

4.5 

1.1 

4.0 

l 

4-178 

1.25 

11.4 

H 

80 

25 

NR 

7.3 

1.7 

5.3 

5-179 

1.26 

12.5 

0 

49 

29 

NR 

6.4 

2.0 

3.8 

Table  5.  White  Sands  Strain  Readings  Rosette  R-2, 
November  28,  1964 


Run 

Mach 

Alt. 
x 1,090 

Vector 

Strain  u In. /In. 
a b c 

Recorded  Pressures  - osf 
M-2  M-7  M-l 6 

2-176 

1.30 

11.9 

A 

NR 

19 

14.3 

3.6 

5.9 

3-177 

1.25 

11.7 

F 

NR 

39 

34 

16.2 

1.1 

4.0 

4-178 

1.25 

11.4 

H 

NR 

52 

53 

19.0 

1.7 

5.3 

5-179 

1.26 

12.5 

D 

NR 

26 

22 

4.8 

2.0 

3.8 

i 

j 
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jki 


North  Wall  East  Wall  South  Wall 


225  -0.707  315  0.707  45  0.707 


t *$&*&  ’^tv^vm  ,*  **t  r it^t  ..y 


P , 

P 

Cl 

North 

P,1 

South 

cos  5 

p.. 

pf 

La. 

pf 

*"10  (“^ 

< 

i fj 

f 

Wall 

Wall 

t 

; 

o, 

1 

6.20000 

3.60000 

3.40600 

•1.00006 

.3663* 

.36636 

-.26061 

-.26091 

i 

2 

6.00000 

3.20000 

4.00000 

••70706 

.66667 

.63633 

-.17606 

-.01666 

3 

3.60000 

2.00000 

3.00000 

•.70700 

.31766 

.53571 

-.26379 

-.27107 

$ 

* 

3.30000 

3.60000 

2.30000 

0.00000 

.66131 

.63366 

-.19260 

•.16255 

1 

3 

3.30000 

6.00000 

2.30000 

0.00000 

.73672 

.63366 

-.12222 

-.36235 

9, 

6 

T. 00000 

3.60000 

2.70000 

0.00000 

.31626 

.36571 

-.2(680 

-.61371 

£ 

T 

3.00000 

3.00000 

0.70000 

•1.00000 

.66667 

1.13536 

-.19107 

•05522 

* 

t 

• 

6.00000 

6.70000 

1.70000 

1.00000 

1.36313 

.13617 

.16663 

• *5079 

y 

Ji 

* 

T. 60000 

6.00000 

1.30000 

.70700 

.66676 

.17165 

-.06630 

.167 

l« 

3.30000 

3.00900 

4.20000 

-.70700 

.36666 

.76263 

-.26713 

-.10.13 

% 

* 

11 

T. 00000 

0.00000 

1.30000 

.70700 

t.2tsn 

.18371 

.10916 

-.71115 

12 

T. 30000 

6.30000 

1.30000 

.70700 

.66661 

.17866 

-.03061 

-.76936 

i 

$ 

13 

6.30000 

3.30000 

1.10000 

•70700 

.61336 

.16623 

-.0(666 

-.77132 

% 

1* 

T. 00000 

6.60000 

3.20000 

0.00000 

.63  M6 

.63716 

-.16216 

-.13995 

I 

U 

3.60000 

6.00000 

2.70000 

0.00000 

.71626 

.66216 

-.1*611 

-.31682 

5; 

{ 

16 

6.20000 

3.60000 

1.S0000 

0.00000 

.36636 

.26166 

-.26091 

-.61610 

3 

IT 

6.20000 

3.60000 

1.30000 

0.00000 

.36636 

.26166 

-.26091 

-.61636 

16 

3.00000 

2.00000 

3.0C000 

-.70706 

.76336 

.67616 

-.12*67 

-.01126 

1* 

3.30000 

6.30000 

1.30090 

.70706 

1.16666 

.26326 

.07506 

-.61031 

20 

3.00000 

2.00000 

4.00001 

-.70700 

.76336 

1.17666 

-.1*667 

.07169 

21 

3.60000 

3.00000 

3.20000 

•l.ooooe 

.33371 

.37163 

-.271(7 

-.26166 

22 

3.30000 

7.00000 

1.30000 

1.00000 

1 .63366 

.26362 

.1363* 

-.36816 

23 

6.20000 

3.20000 

3.70000 

•1.00000 

.31613 

.61613 

-.2672* 

-.61652 

26 

3.10000 

3.60000 

3.00000 

-.70700 

.66667 

.76310 

-.17609 

-.12779 

23 

T, 00000 

3.70000 

1.30000 

.70703 

.72132 

.16636 

-.1*173 

-.7(366 

26 

3.60000 

3.00000 

3.40000 

-.70700 

.33371 

.60716 

-.271(7 

-.21671 

2T 

6.00000 

3.00000 

2.30000 

0.00000 

1.16367 

.66636 

.07123 

-.126*7 

20 

3.20000 

6.00000 

1.20000 

0.00000 

1.13363 

.21077 

.06215 

-.616(2 

20 

1.00000 

3.10000 

3.10000 

-1.00000 

1. 63136 

1.61138 

.21261 

.21261 

30 

6.00000 

1.30000 

2.00000 

0.00000 

.23666 

.33333 

-.602(6 

-.67712 

35 
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respectively.  Columns  6 and  8 contain  the  common  logarithms  of  the 
ratios  from  Columns  5 and  7,  respectively.  It  will  be  shc'/n  later  how 
these  ratios  and  logarithms  are  used  in  the  present  statistical  model. 

3.  EDWARDS  AIR  FORCE  BASE  TESTS 

The  sonic  boom  tests  conducted  at  Edwards  Air  Force  Base. 
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California,  in  June  1966  were  for  measurement  of  both  human  response 
and  structural  reaction.  Like  the  previous  tests,  these  were  conducted 
under  FAA  sponsorship.  Stanford  Research  Institute,  under  the  direc- 
tion of  Karl  Kryter.  performed  the  human  response  stu  ,.e3  a 
John  A.  Blume  and  Associates  again  performed  the  struc  s 

[31, 34]  . NASA  also  participated  in  reduction  and  analysis  o.  ..c  sonic 
boom  data. 

As  part  of  the  Edwards  structural  test  two  test  houses,  shown 
in  Figures  15  and  16,  were  extensively  instrumented.  Of  particular 
interest  in  the  present  study,  five  windows  were  instrumented  with 
strain  gauges  and  two  indoor  microphones  were  installed.  There  was 
also  extensive  instrumentation  of  walls  and  ceilings. 


Between  the  Edwards  tests  and  the  White  Sands  tests,  there  I 

l 

*• 

were  seven  windows  of  various  sizes  instrumented.  These  are  sum-  - 

marized  in  Table  8.  The  strain  reading'?  on  these  windows  form  the 


basis  of  the  present  statistical  model  of  the  dynamic  amplification 

f artnr 
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♦Aluminum  Slidinp  Windows 


See  Fig, 25  for  SG 
locations  on  this  window 


Table  8.  Summary  of  Instrumented  Windows 


Of  the  Edwards  data  there  were  102  data  sets  which  were  found 
to  be  complete  and  suitable  for  use  in  this  study.  These  correspond 
to  102  different  overflights  whose  readings  are  shown  in  Table  9.  Note 
that  for  each  overflight  there  are  strain  gauge  readings  from  windows 
on  five  channels,  internal  pressure  readings  from  two  microphones, 
and  a free  field  pressure  reading.  No  readings  of  external  pressure 
were  made  in  the  Edwards  series.  The  readings  shown  in  Table  9 form 


the  basis  of  the  correlation  analysis  which  led  to  the  statistical  model. 
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Table  9.  Summary  of  Edwards  AFB  Data  (Continued) 


Channel  313  changed  from  colonial  window  to  large  garage  window  from  here  on 


Summary  of  Edwards  AFB  Data  (Continued) 
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IV.  RESPONSE  PROBABILITY 
DENSITY  FUNCTION  TECHNIQUE 

All  the  calculations  of  breakage  probabilities  in  this  report  are 
based  on  the  response  probability  density  function  (pdf)  technique.  In 
order  to  understand  the  analysis  it  is  essential  that  the  concept  of  a 
probability  density  be  fully  understood.  Therefore,  this  chapter  begins 
with  a review  of  that  concept. 

1.  PROBABILITY  DENSITY  FUNCTIONS 

For  illustration,  assume  that  a group  of  402  readings  have  been 
taken,  representing  the  breaking  strengths  of  steel  rods.  The  readings 
are  then  sorted  according  to  size  and  the  number  of  readings  in  each 
250  psi  interval  is  tabulated.  The  number  of  readings  in  each  interval 
is  then  graphed  as  a bar  chart  as  shown  in  Figure  17a.  (If  the  readings 
were  from  glass  rather  than  steel  the  shape  of  the  distribution  would 
be  different,  as  will  be  discussed  later. ) This  particular  form  of 
bar  chart,  showing  the  number  of  readings  in  each  interval,  is  known 
as  a histogram.  The  shape  of  the  histogram  in  Figure  17a  is  the  well- 
known  bell  shape  which  occurs  often  for  distributions  of  men's  heights, 
students'  grades,  and  many  other  measurements. 
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Now  assume  that  for  each  interval  shown  in  Figure  17a  the  num- 
ber of  readings  is  divided  by  the  total  number  of  readings  for  all 
intervals  combined.  This  will  give  the  slightly  different  histogram 
shown  in  Figure  17b.  The  ordinate  01  Figure  17b  shows  the  percentage 
of  the  total  number  of  readings  that  are  in  each  interval.  If  the  sample 
rods  tested  were  typical  of  the  rods  being  manufactured.  Figure  17b 
could  also  be  viewed  as  a histogram  of  probabilities.  Then  it  could 
be  observed  from  the  histogram  that  the  probability  is  .17  of 
picking  a rod  at  random  whose  strength  is  in  the  interval  centered  at 
52.  0 kpsi. 

Note  that  the  histogram  of  Figure  17b  is  a discrete  function. 

For  most  analyses,  it  is  more  convenient  to  work  with  continuous 
functions,  such  as  the  curve  in  Figure  17c.  Note  that  the  curve  looks 
very  much  like  Figure  17b.  In  fact,  it  is  what  would  be  obtained  if 
there  were  a very  large  number  of  samples  and  the  strength  intervals 
were  made  smaller  and  smaller.  The  curve  of  Figure  17c  is  an 
example  of  a probability  density  function.  The  ordinate  of  the  curve 
is  a density  in  probability  per  unit  length.  For  the  normal  curve 
shown  (also  called  a gaussian  density  function),  the  probability 
density  function  is 

1 

f(x)  

y/2tV 


r+m  v~  v <rs  x ,v 


where  n is  the  mean  and  a is  the  standard  deviation.  To  find  the 
probability  that  a stress,  x,  falls  between  two  values,  and  x2. 
merely  take  the  definite  integral  of  the  probability  density  function 
between  x*  and  x2 . Obviously  the  integral  of  a probability  density 
function -from  -®  to  co  is  1 . 

Now  observe  the  probability  density  function  of  a stress  from 
some  form  of  random  loading  as  shown  in  Figure  18.  Note  the  wide 
ranbe  of  values  that  may  be  expected  for  readings  of  stress  from  the 
randomly  occurring  loads  of  this  type.  Assume  that  this  probability 
density  function  has  been  determined  from  a very  large  number  of 
experimental  values.  Now  suppose  that  another  large  sample  of  the 
same  form  of  random  loading  is  taken.  The  probability  density 
functions  should  be  almost  identical  for  the  two  large  samples.  (If 
both  sample  sizes  are  infinite  they  will  be  identical. ) Thus,  although 
individual  stress  values  are  random,  the  shape  of  the  density  function 
is  for  all  practical  purposes  constant.  The  probability  that  stresses 
will  fall  in  a given  range  can  be  predicted  with  certainty. 

Similar  to  the  probability  density  function  for  stress,  a proba- 
bility density  function  for  the  breaking  strengths  of  the  specimens  is 
shown  in  Figure  19,  Again,  the  probability  that  a strength  reading 
w:ll  be  in  a given  range  can  be  determined. 
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Figure  II,  Probability  Ssniity  Function  of  iraking  ftrongth 
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Now,  suppose  one  of  the  stresses  <rmi  from  the  loads  selected 
at  random  acts  on  one  of  the  random  strength  samples  (tq^.  For  this 
case,  an  effective  factor  of  safety  Ne  can  be  defined  by 


N 


°G; 


ek 


(6) 


mi 


Since  the  stress  and  strength  are  assumed  to  be  statistically  indepen- 
dent, the  probability  density  function  of  the  effective  factor  of  safety 
can  be  found  from  the  probability  density  functions  of  the  stress  and 
strength,  The  process  for  finding  the  probability  density  function  of 
Ne  in  the  general  case  is  a convolution -like  procedure  [35]  , which 
can  usually  be  accomplished  by  computer  integration  only, 


F(Ne) 


CD 

=y ^ Ne<rm)fam 

o 


(7) 


Note  that  the  factor  in  the  integral  is  a joint  probability  density  of 
two  variables. 


Considerable  simplification  is  possible  if  the  stress  and  strength 
probability  density  functions  are  both  lognormal.  The  lognormal  prob- 
ability density  function  has  the  form 


f(x) 


xo\/2"ir 


exp 


- -~(ln  x -m)2 
2* 

«a 


(8) 


for  x>0  and  oO,  and  f(x)  = 0 for  x>0.  A variable  is  distributed 
lognormal  if  its  logarithm  is  distributed  normal. 


..m’rsr?? 


If  imposed  stress  and  material  strength  are  both  distributed 
lognormal,  the  probability  density  of  the  effective  factor  of  safety  will 
then  also  be  lognormal.  This  is  because  the  logarithm  of  the  effective 
factor  of  safety  will  be  the  difference  of  the  logarithms  of  the  stress 
and  strength,  and  the  difference  of  two  normally  distributed  random 
variables  is  normally  distributed  with  its  mean  equal  to  the  difference 
of  the  means  and  its  variance  equal  to  the  sum  of  the  variances.  The 
probability  of  breakage  can  then  be  found  as  the  area  under  the  curve 
of  the  probability  density  of  the  effective  factor  of  safety,  for  effective 
factor  of  safety  less  than  1.  Equivalently,  this  is  the  area  under  the 
curve  of  the  density  of  the  logarithm  of  the  factor  of  safety  to  the  left 
of  zero. 

The  probability  density  functions  for  the  effective  factor  of  safety 
and  its  logarithm  are  shown  in  Figure  20.  The  easiest  way  to  find  the 
probability  of  breakage  in  this  case  is  to  work  with  the  normal  curve 
of  the  logarithm  of  Ne  in  Figure  20b,  since  the  probabilities  under  a 
normal  curve  are  a tabulated  function. 

The  following  sections  will  develop  the  probability  density  func- 
tions of  sonic  boom  stress  and  material  strength  so  that  the  above 
nique  can  be  utilized. 
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A.  PROBABILITY  DENSITY  FUNCTION  OF  EFFECTIVE  FACTOR 


LOGARITHM  OF  EFFECTIVE  FACTOR  OF  SAFETY 


B.  PROBABILITY  DENSITY  FUNCTION  OF  THE 

LOGARITHM  OF  THE  EFFECTIVE  FACTOR  OF  SAFETY 


Figure  20.  Probability  Density  Function  of  the  Effective  Factor 
of  Safety  for  the  Lognormal  Case 


53 


■£  fit  wx  gfc*  M&Vri  ft*  &!&&*****&*!£ 


:r^as»Hs 


2.  SONIC  BOOM  STRESS 

With  the  definit'ons  of  the  peak  overpressures  and  the  DAF 
discussed  in  Chapter  II,  a correlation  analysis  was  performed  using 
experimental  data  from  the  sonic  boom  tests  performed  at.  White  Sands 
in  1965  [29,30,33]  and  Edwards  Air  Force  Base  in  1966  [31 , 32,  33  ] . 
These  data  consisted  of  values  of  pf,  pe,  and  p^  measured  from  micro- 
phones and  concomitant  values  of  um  measured  with  strain  gauges  on 
the  windows  on  the  test  houses. 

Correlation  analysis  [1.4]  of  the  Edwards  data  indicated  that 
strain  in  the  test  windows  correlated  well  with  free  field  pressure  and 
poorly  with  other  variables.  The  correlation  coefficients  calculated 
between  pf  and  the  strains  on  four  windows  were  0.66,  0.  69,  0.  69, 
and  0.42.  All  these  results  indicate  significant  correlation  at  the 
99  percent  confidence  level.  Correlation  coefficients  were  also  com- 
puted between  the  strains  and  the  impulse  PfT/2,  but  correlations 
were  found  to  be  less  than  for  the  pressure  pf.  Similarly  correlations 
with  wave  energy  pj  T/3  were  also  found  to  be  less  than  those  with  p^. 
These  results  indicate  that  the  free  field  pressure,  pf,  is  the  correct 
variable  to  use  in  formulating  a linear  expression  for  stress. 

Three  variables  that  were  discarded  on  the  basis  of  low  correla- 
tion were  the  internal  pressure  pjf  the  rise  time  of  the  free  field 
pressure  waveform  tr,  and  the  sonic  boom  duration  T.  The  strains 
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on  two  windows  were  found  to  have  correlations  of  -0.  06  and  +0. 18 
respectively  with  Pj/pf.  These  low  correlations  with  signs  in  both 
directions  suggest  that  internal  pressure  is  for  ail  practical  purposes 
uncorrelated  with  strain.  The  DAF's  on  five  windows  were  found  to 
have  correlations  with  rise  time  of  +0.20,  -0.28,  +0.01,  and  -0.06. 
This  indicate."  that  DAF  is  not  correlated  with  rise  time.  Correlation 
analysis  and  examination  of  scatter  diagrams  also  showed  that  stress 
is  unrelated  to  the  sonic  boom  duration  T except  in  the  case  where  a 
very  short  duration  boom  impinges  on  a very  large  storefront  window. 
In  this  exceptional  case,  such  as  might  occur  from  an  F-104  boom, 
there  is  a lowpass  filter  effect.  The  short  duration  boom  is  ineffec- 
tive in  exciting  the  window  above  its  natural  frequency  and  a low  DAF 
occurs,  perhaps  of  the  order  of  0.  5. 

A correlation  of  0.78  was  found  between  the  peak  external 
pressure  pe  and  the  peak  free  field  pressure  Of.  Scatter  diagrams 
of  pe/Pf  from  the  White  Sands  data  indicated  a marked  dependence  on 
the  cosine  of  the  flight  path  angle,  0.  On  this  basis  a regression 
analysis  was  performed  resulting  in  the  relation 

log10(pe/pf)  = .1427  cos  9 - .1258.  (S) 

Having  observed  the  results  of  the  correlation  analysis  it  was 
now  desired  to  formulate  a statistical  rel.  .ionship  between  the  stress 
and  the  nominal  overpressure.  Data  from  the  Oklahoma  City  sonic 
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boom  tests  of  1964  [28]  indicated  a correlation  of  0.42  between  the 
nominal  pressure  pQ  and  the  free  field  pressure  p^.  The  dependence 
of  external  pressure  pe  and  free  field  pressure  pf  had  also  been  noted 
earlier  as  had  the  dependence  of  stress  on  pf.  Thus  the  most  suitable 
form  for  the  stress  relation  was  determined  to  be 

<7m  = PoWPo^Pe/Pf^d/Pe^m^d*-  (10) 

In  the  above  form  all  the  ratios  shown  are  measurable  values  for 
which  data  are  available.  The  Oklahoma  City  tests  provided  3500  values 
of  Pf/pQ  and  the  White  Sands  tests  provided  900  values  of  Pe/Pf-  The 
value  o-^/pg  is  deterministic;  it  can  be  visualized  as  a stress  factor  F 
such  as  that  given  in  Equation  3 for  the  flat  plate  case.  The  value 
ffm/^d  is  the  D^F>  f°r  which  500  values  are  available  from  the 
Edwards  Air  Force  Base  and  White  Sands  tests. 

It  is  readily  apparent  from  the  physical  situation  that  the  ran- 
dom variable  factors  in  Equation  10  are  ail  statistically  independent. 

The  value  of  (pf/pD)  depends  on  the  inhomogeneities  in  the  atmos- 
and  on  the  deviation  of  the  aircraft  from  its  planned  altitude  and  Mach 
number.  The  random  variable  (pe/pf)  depends  on  the  geometry  of  the 
surfaces  involved  and  their  reflection  coefficients.  The  DAF  depends 
on  window  size  and  mounting.  Since  the  three  random  variables 
depend  on  independent  phenomena,  they  can  be  expected  to  be 


independent. 


Using  the  data  from  the  Oklahoma  City,  White  Sands,  and 
Edwards  Air  Force  Base  tests,  histograms  were  computer  plotted 
for  Pf/pQ.  Pe/Pf-  an<*  OAF.  The  results  [1,4]  showed  that  their 
pdf's  were  to  a reasonable  approximation  lognormal.  Thus  the  histo- 
grams of  the  logarithms  appeared  normally  distributed  as  did  the 
example  of  (Pf/pQ)  oiiown  m Figure  21;  a gaussian  curve  of  the 
same  mean  and  variance  is  also  plotted  for  comparison.  Thus  we 
have  verified  that  the  stress  was  indeed  lognormal.  The  next 
section  considers  the  probability  density  function  of  the  material 
strength. 

3.  MATERIAL  STRENGTH 

The  breaking  stress  cq  of  a materiel  can  be  expressed  in  terms 
of  a breaking  pressure  Pq  multiplied  by  a stress  factor  F,  which 
depends  on  the  loading  configuration. 

°G  = Pgf-  (11) 

The  materials  we  are  discussing  here  are  brittle  in  nature,  such  as 
glass,  plaster,  and  mortar.  For  each  of  these  materials  the  ques- 
tion arises  as  to  what  the  shape  of  the  pdf  of  the  breaking  pressure 
is. 

We  have  been  able  to  obtain  data  on  the  shape  of  the  breaking 
pressure  for  some  of  the  materials  and  found  them  to  be  lognormal. 


THE  VALUE  GIVEN  fO»  tHE  INTERVAL  »ioth  is 
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Figure  21.  Experimental  Probability  Density  Function  of  Logarithm  of 
External  Pressure /Nominal  Pressure,  Station  3 


In  our  previous  work  on  glass  [1,4]  we  obtained  a large  data  base 
of  glass  breaking  pressure  readings  from  the  Libbey -Owens -Ford 
Company  [36]  . By  plotting  the  histograms  of  the  breaking  pressures 
we  found  the  probability  density  function  of  the  glass  breaking  pres- 
sure to  indeed  be  lognormal.  Similarly  we  obtained  experimental 
data  on  the  tensile  strength  of  mortar  [37,  38]  and  determined  that 
the  pdf  is  also  lognormal.  For  other  materials  we  have  been  unable 
to  find  sufficient  experimental  strength  data  to  plot  a histogram  indica- 
tive of  a pdf  but  on  the  basis  of  their  brittle  nature  and  their  similarity 
to  glass  or  mortar  we  have  assumed  their  pdf's  to  be  lognormal 
also. 

4.  PROBABILITY  OF  BREAKAGE 

Having  examined  the  pdf's  of  the  stress  and  the  strength  and 
found  them  to  be  lognormal,  the  response  probability  density  function 
technique  can  now  be  applied  to  find  the  probability  of  breakage.  The 
effective  factor  of  safety  Ng  is  given  by 

Ne  = aC>m-  (12) 

where  <Tq  is  the  breaking  pressure  of  the  material  and  am  is  the  maxi- 
mum stress  created  by  the  sonic  boom.  Using  Equations  10  and  11 
for  the  maximum  sonic  boom  stress  orm  and  the  material  strength 
<7q  we  have  after  taking  logarithms 
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loSlO  ffm  = loeiO  Po  + lQglO  <Pf/Po>  + l°giC  (Pe/pf)  <13> 

4 • 210®*  + logio(E)AF), 

logio  o'G  = logl0F  + log10  Pg-  (14) 

Thus  there  results  from  Equations  13  and  14 

E Jlogio  Ne|  = E jlogio  Pg|  ~ lQglO  Po  " E jl°gl0  (PfW|  (15> 
- E ]log10(pe/pf)j  - E jlogi0(DAF)j. 

and 

Var  (log10Ne)  = Var  jlogio  (pf/p0)J  + Var  [log10(pe/pf)J  (16) 
+ Var  ^logio(DAF)j  + Var  jlogio  (PG>j 
The  pdf's  of  the  terms  involved  are  illustrated  in  Figure  22. 


The  mean  and  variance  of  the  distribution  of  each  logarithm 
were  found  for  each  case.  Note  that  the  mean  must  be  found  after 


taking  logarithms,  since  the  mean  of  the  log  is  not  equal  to  the  log  of 
the  mean  [39]  . For  logjQ  (pf/p0)  the  mean  was  found  to  be  0.  0471 
and  the  variance  was  0,  0446.  For  sonic  booms  equally  likely  from 
any  flight  path  direction*  the  mean  of  logiQ  (pe/pf)  was  -0.1251  and 
its  variance  was  0.0439.  If  instead  only  head-on  booms  are  consid- 
ered then  the  mean  of  log^Q  (Pe/Pf)  is  0.  0174  and  its  variance  is 
0.0309.  For  the  remaining  random  variables  DAF  and  Pq  the  mean 
and  variance  depend  on  the  material  configuration.  The  other  two 
numbers  F and  pQ  are  deterministic. 
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Development  of  Probability  Density  Function  of  log^g 


E 


Having  found  the  mean  and  variance  of  log^o  Ne,  the  probability 
of  breakage  follows  directly.  By  finding  the  value  of 

z = E Jlog10  NeJ  I [Var  (log10  Ne)]y’  (17) 

the  area  to  the  left  of  zero  in  Figure  20  can  be  found  by  entering  a 
standard  gaussian  table. 

This  is  the  technique  which  will  be  used  in  succeeding  chapters 
to  find  the  probability  of  breakage  for  various  materials.  The  next 
chapter  will  consider  glass,  and  subsequent  chapters  will  analyze 
plaster,  bric-a-brac,  and  brick. 


V.  GLASS 


Unlike  metals,  the  microstructure  of  glass  is  amorphous 
rather  than  crystalline.  Because  of  this  fact,  the  practical  design 
strength  of  glass  is  a surface  condition  property  rather  than  a con- 
stant material  property.  For  example,  glass  that  is  acid -etched  and 
coated  with  lacquer  to  protect  its  surfaces  can  have  an  average 
strength  of  250  kpsi.  If  the  same  glass  is  severely  sandblasted, 
its  average  strength  is  reduced  to  2 kpsi  [40]  . Because  of  this 
dependence  on  surface  scratch  condition,  strength  tests  on  seemingly 
identical  panes  of  glass  will  show  a wide  scatter  of  values.  Even  if 
the  panes  of  glass  have  received  the  same  handling,  the  depth  and 
location  of  minute  surface  scratches  will  be  random.  The  exact 
valu*-  ol  breaking  strength  for  a given  pane  cannot  be  predicted; 
however,  the  mean  and  variance  of  large  lots  of  glass  of  the  same 
size,  type,  and  surface  condition  will  be  repeatable. 

In  addition  to  the  variations  due  to  surface  scratch  condition, 
there  are  also  variations  with  loading  geometry,  loading  rate,  ?tmos- 
pheric  moisture  content,  and  composition.  Glass  exhibits  a property 
known  as  "static  fatigue"  in  that  it  is  weaker  for  loads  of  longer  dur- 
ation. Static  fatigue  is  a stress  corrosion  type  of  process  that  takes 
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place  when  water  vapor  gets  into  cracks  during  loading  and  aids  their 
growth  [41]  . Thus  for  sonic  boom  loading,  which  has  a duration  of 
the  order  of  0, 1 sec,  the  strengths  of  glass  will  be  roughly  twice 
those  obtained  in  typical  laboratory  measurements  whose  loading 
duration  is  of  the  order  of  60  sec,  as  shown  in  Figure  23. 

In  our  study  we  were  fortunate  in  being  able  to  use  a data  base 
of  unpublished  static  test  results  on  2500  lites  of  glass,  which  the 
Libbey-Owens-Ford  Company  graciously  made  available.  The  data 
consisted  of  119  sets  representing  a full  range  of  window  sizes  with 
approximately  25  readings  of  breaking  pressure  per  set.  Each 
reading  was  obtained  by  loading  the  glass  incrementally  to  failure  on 
a vacuum  test  fixture,  with  the  effective  load  duration  being  approxi- 
mately one  minute. 

Two  adjustments  were  necessary  in  extrapolating  from  the 
Libbey-Owens-Ford  data  of  new  glass  under  static  loading  to  the 
case  of  used  glass  in  people's  homes  under  sonic  boom  loading. 

First,  the  glass  must  be  considered  to  have  twice  the  strength  under 
sonic  boom  loading  as  under  60  second  laboratory  leading  because  of 
the  static  fatigue  effect.  Counterbalancing  this,  a reduction  of  the 
strength  must  be  allowed  for  the  deterioration  of  the  surface  condi- 
tion over  time.  After  the  glass  emerges  from  the  factory  it  suffers 
scratches  in  shipping,  unpacking,  handling,  installing,  cleaning,  and 
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weathering.  Although  the  glass  is  weakened  by  this  progressive 
scratching  over  time,  there  appears  to  be  no  cumulative  dynamic 
fatigue  effect  due  to  repeated  sonic  boom  loading  [42]  . On  the  basis 
of  comparing  the  only  available  published  strength  data  on  used  glass 
with  test  results  on  the  corresponding  size  of  new  glass  [42]  , it  was 
observed  that  the  used  glass  appeared  decreased  in  strength  by  a 
factor  of  2,  as  shown  in  Figure  24.  Thus  the  two  factors  of  2 above 
compensate  ror  each  other,  and  the  strength  of  new  glass  under  static 
loading  can  be  assumed  to  be  the  same  as  that  of  used  glass  under 
sonic  boom  loading. 

Computer-generated  histograms  [1,4]  were  run  for  the  glass 
breaking  strength  and  the  glass  breaking  pressure.  Both  were  found 
to  be  lognormal  to  a reasonable  approximation. 

It  was  necessary  to  provide  in  the  statistical  model  for  the 
presence  of  already  cracked  glass  in  the  window  population.  Although 
there  appear  to  be  no  published  test  data  on  cracked  glass,  experi- 
menters in  the  glass  community  have  said  that  10  percent  of  the 
strength  of  ordinary  used  glass  is  a good  rule  of  thumb.  On  this 
basis  the  present  study  has  assumed  that  the  mean  strength  of 


cracked  glass  is  one-tenth  that  of  used  glass,  and  that  the  variance 


of  the  distribution  of  the  logarithm  of  its  strength  is  the  *,  >>  as  that  | 

3 

of  used  gla^j-.  A previous  study  [31  ] found  cracked  glass  to  be  ' i 


.61  percent  of  the  population.  On  this  basis  the  glass  population 


was  assumed  in  this  study  to  consist  of  89.39  percent  gooa  glass  and 
0.61  percent  cracned  glass.  The  pdf  of  the  logarithm  of  strength 
for  the  entire  population  is  shown  synthesized  in  Figure  25.  Note 
that  in  his  model  all  glass  is  assumed  to  be  either  used  glass  in 
good  condition  or  cracked  glass.  No  provision  has  been  made  for 
brand  new  glass  or  for  any  state  between  cracked  and  good.  Using 
this  model  the  probabilities  of  breakage  for  the  healthy  glass  Ph  and 
the  cracked  glass  P<y  can  then  be  combined  into  a population  proba- 
bility of  breakage  P by  using  the  relation 

P = 0.9939PH  + 0.0061  PC.  (18) 


In  applying  the  response  pdf  technique  To  predicting  the  sonic 
boom  breakage  probability  in  a group  of  windows,  all  the  windows 
will  have  the  same  pdf's  of  logig  (Pf/p0)  anc*  -°gl0  (pe/PfB  But  the 
pdf's  of  logjQ  iDAF)  and  Iogjg  ((Vj)  wtli.  in  general,  be  different  for 
eav.n  window  depending  on  its  geometry  . Sev  m ivpical  window  types 
for  which  data  were  available  are  indicated  {ri  Table  8.  In  t .is 
table  the  column  headed  fQ  gives  the  natural  frequency  of  the  tite, 
which  is  important  only  in  the  case  where  a short  duration  boom 
excites  a very  large  window.  The  sixth  column  gives  the  stress 
factor  F trom  Equation  3.  The  last  column  gives  the  value  of 
(abj'^/h;  stress  is  roughly  p.oporlioral  to  the  square  of  this 


Model  of  the  Probability  Density  Function  of  the  Logarithm  of 
the  Breaking  Strength  of  the  Glass  Population 


parameter  for  small  deflections.  The  means  and  variances  of  the 
various  factors  for  these  windows  for  all  flight  paths  combined  and 
for  head-on  booms  are  shown  in  Tables  10-13. 

Note  that  for  a given  window  size  the  mean  of  logigNg  °*'  Equa- 
tion 15  decreases  with  increasing  pQ.  Probability  of  breakage  vs. 
nominal  overpressure  p0  was  plotted  for  the  seven  typical  window 
sizes  in  good  condition;  the  results  were  straight  lines  on  lognormal 
paper.  These  are  shown  f ir  the  cases  of  booms  equally  likely  from 
any  direction  and  for  head-on  booms  in  Figures  26  and  27  respectively. 
These  results  agree  well  with  sonic  boom  damage  claims  data  at  low 
overpressures  of  0.6  breaks  per  million  panes  boomed  at  a nominal 
overpressure  of  1 psf  [43. 44]  . The  results  also  agree  with  NASA 
experiments  at  high  overpressures  [45]  . 

Probabilities  of  breakage  were  also  computed  for  a model  win- 
dow population  consisting  of  both  good  glass  and  cracked  glass  using 
Equation  17.  The  results  are  shown  in  Tables  14-21  for  various 
nominal  overpressures  from  1 to  100  psf  with  booms  equally  likely 
from  any  direction  or  head-on.  The  composite  glass  population 
probabilities  assume  equal  numbers  of  each  of  the  seven  window  types. 
Note  that  for  this  low  orpressure  the  majority  oi  breaks  are  likely 
to  occur  in  already  cracked  windows. 
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Lble  10.  Mean  of  the  Logarithm  of  the  Efi'3ctive  Factor 
of  Safety  for  All  Flight  Paths  Combined 


Table  12.  Means  of  the  Logarithm  of  the  Effective  Factor 
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Table  14.  Probability  of  Glass  Breakage  Calculated  for  Model  Window  Population  for 
Ail  Flight  Paths  Equally  Likely  at  a Nominal  Overpressure  of  1 psf 


Window  Type 


Glass  Population  P = 67.6X10' 
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Glass  in  Good  Condition 

Cracked  Glass 

Combined 

Window 

No. 

Mean 

lo*10Ne 

z 

P 

Mean 

logioN< 

Z 

a 

• 

P 

P 

1 

1.8836 

5.  67 

.7  X10"4 

. 8836 

2.  65 

0040 

2.  44X10"5 

2 

1.9250 

5.71 

. 6 X10-1 

. 9250 

2.  74 

.0031 

1.  89X10"’ 

3 

1.8892 

5.  59 

1.  1 X10"‘ 

. 8892 

2.  63 

.0043 

2.  62X10"* 

4 

1. 6875 

4.  47 

4.  3 X10"4 

. 6875 

1.82 

. 0346 

21.  53X10"* 

5 

2.0821 

5.  11 

1.  67X10"’ 

1.0821 

2.  65 

. 0040 

2.46X10"* 

6 

1.7637 

4.  51 

3.4  X10"‘ 

.7637 

1.  95 

. 0254 

15.  83X10”* 

H 

1 

2.3107 

5.41 

. 3 X10"’ 

1.3107 

3.07 

.0010 

. 61X10“* 

Good  Glass 

V 

1.  1X10"4 

Cracked  Glass 

tl 

&P 

0.0109 

jag  os  t aaa  ■ •- — 


Table  15.  Probability  of  Glass  Breakage  Calculated  for  Model  Window  Population 
for  Head-On  Flight  Path  at  a Nominal  Overpressure  JlZ P°PUlatl°n 


Window 

No. 


Glass  in  Good  Condition 
Mean 

1O*10Ne  " P 


Cracked  Glass 
Mean 

loho%  2 


Window  Type 
Combined 


1 

1.7411 

5.  59 

1.  2X10“* 

.7411 

2.  38 

. 0087 

5.13X10' 

2 

1.7825 

5.  62 

9.  8X10"* 

. 7825 

2.47 

.0068 

4.  15X10“J 

3 

1. 7467 

5.49 

2.  1X10“* 

.7467 

2.  34 

. 0096 

5.  86X10"* 

4 

1.5450 

4.  29 

8.  6X10"* 

. 5450 

1.  51 

. 0655 

40.  81X10' 5 

5 

1.9396 

4.  96 

3.  5X10’’ 

.9396 

2.40 

. 0082 

5.  04X10“ 5 

6 

1.  6212 

4.  34 

6.  9X10“* 

.6212 

1.  66 

. 0485 

30.  27X10~S 

7 

2. 1682 

5.  26 

7.4X10"*  1 

. 1682 

2.  83 

. 0023 

1.  41X10-5 

Good  Glass 

PH 

= 2. 

3X10~‘ 

Cracked  Glass 

pc 

= .0214 

Glass  Population  p 

= 1. 

33X10'4 
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Table  16.  Probability  of  Glass  Breakage  Calculated  for  Model  Window  Population  for 
All  Flight  Paths  Equally  Likely  at  a Nominal  Overpressure  of  2 psf 


Window  Type 

Glass  in  Good  Condition  Cracked  Glass  Combined 


Window 

No. 

Mean 

log,  .N 
610  e 

z 

P 

Mean 

lo*10N'e 

z 

P 

P 

1 

1.5826 

4.76 

1.02X10“* 

.5826 

1.  75 

. 0401 

2.46X10"* 

2 

1. 6240 

4.  82 

. 80X10"* 

. 6240 

1.  85 

. 0322 

1.  97X1 0"4 

3 

1.5882 

4.70 

!.  31X10"* 

. 5882 

1.  74 

. 0409 

2.50X10"* 

4 

1.3865 

3.  67 

1.  3 X10"* 

.3862 

1.  02 

. 1539 

1.07X10"3 

5 

1.7811 

4.  37 

5.  5 X10“* 

.7811 

1.  92 

.0274 

1.  73X10"*. 

6 

1.4627 

3.  74 

1.  0 X10“* 

.4627 

1.  w 

. 1379 

9.41X10"* 

7 

2.0097 

4.  70 

1.  31X10"* 

1.0097 

2.  36 

. 0081 

5.  07X10"’ 

Good  Glass 

PH 

= 3. 

43X10“* 

Cracked  Gla 

ss  Pc 

* 0. 

0629 

Glass  Population  P = 4.  18X10“* 
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i AilJ±>£jLfA  **«<*i^Mt* 


1 

1.4401 

4.  63 

2.  0 X10"‘ 

. 4401 

1.41 

. 0793 

4.  86X10'4 

2 

1.4815 

4.  70 

1.4  X10-4 

. 4815 

1.  52 

.0643 

3.  94X1  j‘* 

3 

1.  4457 

4.  54 

2.  8 X10~* 

.4457 

1.40 

. 0808 

4.  96X1 0~* 

4 

1.2440 

3.  45 

2.  81X10-4 

. 2440 

. 68 

. 2451 

1.77X10'3 

5 

1.6386 

4.  13 

1.5  X10‘* 

.6386 

1.63 

. 0515 

3.  29X10'4 

6 

1.3202 

3.53 

2.20X10'4 

.3202 

.86 

. 1949 

1.  41X10~J 

7 

1.8672 

4.53 

2.8  X10"‘ 

.8672 

2.  10 

. 0179 

1.  12X10  4 

Good  Glass  PH  = 7.  5 x 10~5 

Cracked  Glass  Pc  = ,105 

Glass  Population  P = 7.15  x 10~4 


* f »)^f  <%» 


Table  18.  Probability  of  Glass  Breakage  Calculated  for  Model  Window  Population 
for  Head-On  Flight  Path  at  a Nominal  Overpressure  of  4 psf 

Window  Type 


Window 

No. 

Glass  in 

Good  Condition 

Cra 

eked  Glass 

Combined 

Mean 

l°SioNe 

z 

P 

Mean 

l°8lONe 

z 

P 

P 

1 

1.  1390 

3.  66 

1.  4X10"4 

. 1390 

.45 

. 3264 

2.  13X10'1 

2 

1. 1804 

3.  72 

. 9X10"* 

. 1804 

.57 

. 2843 

1.  82X10'1 

3 

1.  1446 

3.  60 

1.  6X10”4 

. i446 

.45 

. 3264 

2.  15X10* 3 

4 

. 9430 

2.  62 

4.4X10"3 

-.0570  - 

16 

. 5836 

7.  81X10-3 

5 

1.  3375 

3.  42 

3.  2x:t4 

. 3375 

. 86 

. 1949 

1. 50X10-3 

6 

1.0191 

2.73 

3.  2X10 

.0191 

. 05 

.4801 

6.  11X10"3 

7 

1.5661 

3.  80 

. 8X10"" 

.5661  1 

. 37 

. 0853 

. 60X1  O'3 

Good  Glass  P « 1.20X10 

ri 

Cracked  Glass  P^,  - .323 

Glass  Population  P = 3.  16X10 
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Table  19.  Probability  o i Glass  Breakage  Calculated  for  Model  Window  Population 
for  Head-On  Flight  Path  at  a Nominal  Overpressure  of  20  psf 


Window  Type 


Window 

No. 

Glass  in 

Good  Condition 

Cracked  Glass 

Combined 

Mean 

loglONe 

2 

P 

Mean 

lc>gl0NTe 

z 

P 

P 

1 

.4401 

1.41 

.0793 

-.5599 

-1.  80 

. 964. 

. 085 

2 

.4815 

1.52 

.0643 

-.5185 

-1.63 

. 9485 

.070 

3 

.4457 

1.40 

.0808 

-.5543 

-1.  74 

. 9591 

. 086 

4 

.2440 

.68 

. 2451 

-.7560 

o 

i 

. 9821 

.250 

5 

.6386 

1.63 

.0515 

-.3614 

- .92 

. 8212 

. 056 

6 

.3202 

.86 

. 1949 

-.6798 

-1.82 

.9656 

. 200 

7 

.8672 

2.  10 

.0179 

-. 1328 

- .32 

.6255 

. 022 

Good  Glass 

PH  ' 

. 105 

Cracked  Glass 

pc  = 

. 895 

Glass  Population 

P = 

. 110 
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-3 


1 

. 1-JO 

.45 

.3264 

• .8610 

-2.76 

.9971 

.330 

2 

.1804 

.57 

.2843 

• .8196 

-2.58 

.9951 

.289 

3 

.144* 

.45 

.3264 

- .8554 

-2.69 

.9964 

.330 

4 

-.0570 

- .16 

.5636 

-1.0570 

-2.93 

.9983 

.566 

5 

.3375 

.86 

. 1949 

- .6625 

-1.69 

.9545 

.200 

6 

.0191 

.05 

.4801 

- .9809 

-2.62 

.9956 

.483 

7 

.5661 

1.37 

.0853 

- .4339 

-1.05 

.8531 

.090 

Good  Glass 

PH  * 

.323 

Cracked  Glass 

II 

u 

(U 

.970 

Glass  Population 

p = 

.327 

Table  21.  Probability  of  Glass  Breakage  Calculated  for  Model  Window  Population 
for  Head-On  Flight  Path  at  a Nominal  Overpressure  of  100  psf 


Glass  in  Good  Condition 


Window  Type 

Cracked  Glass  Combined 


Window 

Mean 

Mean 

No. 

IO«10Ne 

z 

p 

1Og10Ne 

P 


1 

-.2589 

- .83 

.7967 

-1.2589 

-4.04 

. 9999 

.798 

2 

-.2175 

- .69 

.7549 

-1.2175 

-3.83 

.9999 

.756 

3 

-.2533 

- .80 

.7881 

-1.2533 

-3.93 

. 9999 

.789 

4 

-.4550 

-1.26 

.8962 

-1.4550 

-4.  04 

. 9999 

.897 

5 

-.0604 

- .15 

.5596 

-1.0604 

-2.71 

.9966 

.562 

6 

-.3788 

-1.01 

.8438 

-1.3788 

-3.68 

.9999 

.845 

7 

+. 1682 

+ .41 

.3409 

- .8318 

-2.02 

.9783 

.345 

Good  Glass 

PH  = 

.711 

Cracked  Glass  P^  = . 9963 
Glass  Population  P = .713 
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The  external  overpressures  are  somewhat  higher  for  head-on 
booms  than  for  flight  paths  equally  likely  from  any  direction  as 
indicated  in  Equation  9.  For  a 1 psf  head-on  boom,  the  breakage 
probability  calculated  for  the  model  glass  population  was  133  breaks 
per  million  panes  boomed;  for  2 psf  head  on,  715  breaks  per  million. 
Note  that  these  values  of  P are  almost  twice  those  shown  in  Tables  14 
and  16  for  the  same  overpressures  and  flight  paths  equally  likely 
from  any  direction.  Breakage  probabilities  for  the  model  population 
were  also  computed  for  head-on  booms  at  higher  overpressures;  at 
4 psf,  3,160  breaks  per  million;  at  20  psf.  110,  000  breaks  per 
million;  at  40  psf.  327,000  breaks  per  million;  and  at  100  psf. 

713,000  breaks  per  million. 

This  concludes  the  analysis  of  glass;  the  remaining  chapters 
extend  the  response  probability  density  function  technique  to  plaster. 


bric-a-brac  and  brick. 
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VI.  PLASTER 


Plaster,  like  glass,  is  a brittle  material  which  is  sometimes 
subject  to  sonic  boom  damage.  Plaster  is  manufactured  by  heating 
(calcining)  gypsum  at  300°-350°F.  The  calcining  process  drives  off 
water  vapor  and  changes  the  state  of  the  material  from  dihydrate 
calcium  sulfate  to  hemihydrate.  The  calcined  gypsum  can  be  used 
to  form  various  plasters  at  the  building  site,  depending  on  the  aggre- 
gate with  which  it  is  mixed.  In  older  construction  the  aggregate  con- 
sisted of  lime  and  sand.  Other  possible  aggregates  are  wood  fiber, 
vermiculite,  and  perlite.  When  the  plaster  and  aggregate  are  mixed 
with  water  they  form  a slurry  and  entrapped  air  bubbles  float  out. 

As  the  plaster  sets,  a crystallization  process  takes  place  and  the 
gypsum  returns  to  its  dihydrate  form  and  bonds  in  the  aggregate 
materials. 


The  proportioning  of  ingredients,  the  thoroughness  of  mixing, 
and  the  removal  of  air  bubbles  ail  depend  on  the  workmanship  of  the 
individual  plasterers.  For  this  reason  plaster  shows  considerable 
variation  in  strength. 


The  strength  of  the  plaster  will  depend  to  a large  extent  on 
whether  the  correct  amount  of  water  is  mixed  in;  too  much  or  too 
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little  will  lower  the  strength.  The  key  parameter  in  determining  the 
mechanical  properties  of  plaster  is  the  water-to-piaster  ratio,  which 
is  called  the  consistency.  If  too  much  water  is  added  the  excess  will 
run  off,  leaving  the  set  plaster  weak.  Iv  too  little  is  added  incomplete 
hydration  will  take  place  and  the  resulting  plaster  will  have  places 
that  have  never  gone  into  the  crystalline  form.  Figures  28  and  29 
[46,48,49]  show  the  variation  in  mechanical  properties. 


There  is  also  a substantial  variation  in  tensile  strength  of 
plaster  depending  on  the  proportion  of  aggregate  and  the  type  of 
aggregate  used.  This  variation  is  indicated  in  Table  22  [50].  Note 
from  the  table  that  plaster  is  considerably  stronger  in  compression 
than  in  tension,  so  that  for  sonic  boom  push-pull  loading  only  tension 
failures  would  be  expected.  From  Table  22  it  appears  that  a typical 
conservative  figure  to  use  for  the  tensile  strength  of  plaster  is 
100  psi.  However,  the  strength  could  get  up  to  350  psi. 


Usually  plaster  is  applied  to  a wall  or  ceiling  in  three  coats. 
The  first  coat  or  "scratch  coat"  is  applied  directly  to  the  lathing. 

Its  surface  is  then  scratched  with  a rough  tool  to  provide  better 
adherence  of  the  following  coat.  The  second  coat,  called  the  "brown 
coat,  " is  somewhat  thicker.  The  "finish  coat"  which  goes  over  the 
brown  coat  is  the  one  that  is  visible  to  the  eye.  It  is  often  composed 
of  "gauging  plaster,  " a combination  of  calcined  gypsum  and  lime. 
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Table  22.  Strength  of  Basecoat  Plaster  for  Mixtures 
With  Various  Aggregates 


Type  of  Plaster 
and  Miv* 

Compressive  Strength 
Range,  Dry,  psi 

Tensile  Strength 
Range,  psi 

Mill'>  <?ed  Perlite 

Piaster' 

600-800 

125-160 

Wood  Fiber  Plaster 

Neat 

150C  2000 

320-330 

Sanded  100:1 

1200-1600 

180-240 

O/psum  Neat  Piaster 

V.  <e<*  With: 

Sw  i00:2 

750-1100 

140-170 

100:2-1/2 

650-850 

120-145 

100:3 

550-750 

100-125 

Perlite  100:2 

600-800 

120-150 

100:3 

450-600 

95-110 

Vermiculite  100:2 

400-525 

90-100 

100:3 

250-325 

75-85 

^Aggregates  shown  in  cu  ft,  piasters  in  lb 

Note:  Average  laboratory  test  results;  actual  job  strengths  may 
vary  from  these  data  because  of  job  conditions  and 
methods  of  mixing  and  handling. 


The  total  thickness  for  all  three  coats  combined  is  usually  1 /2 


to  3/4  inch. 


The  analysis  of  plaster  failures  on  buildings  is  much  more 


difficult  than  those  of  glass,  because  plaster  never  appears  as  a lone 


element  subjected  to  sonic  boom  as  glass  does.  Plaster  is  always 


used  in  conjunction  with  a wall  or  ceiling  configuration,  which  sup- 


plies the  vast  majority  of  the  load-resisting  capacity  of  the  assembly. 


Each  wall  or  ceiling  configuration  thus  represents  a different  support 


for  the  plaster  mounted  upon  it  and  the  stress  on  the  plaster  can  be 


expected  to  vary  accordingly.  The  more  the  support  structure  de- 


flects, the  greater  the  stress  on  the  plaster.  Because  of  the  many 


variations  possible,  the  best  that  can  be  attempted  within  the  scope 


of  this  study  is  to  perform  the  analysis  for  a few  configurations  which 


seem  typical. 


Besides  the  structural  variation  another  problem  which  compli- 


cates piaster  failure  predictions  is  the  nature  of  plaster  failures. 


Generally  the  overstressing  of  building  plaster  in  tension  merely 


results  in  hairline  cracks,  but  not  in  immediate  catastrophic  failure. 


This  is  much  different  from  sonic  boom  failure  of  glass  where  it  is 


like;y  that  the  window  will  shatter,  or  that  a highly  visible  crack 


will  appear.  Plaster  cracks  by  comparison  are  very  thin  and  are 
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difficult  for  the  untrained  observer  to  detect  with  the  naked  eye, 
especially  if  they  occur  in  a ceiling.  Only  when  a catastrophic  failure 
occurs,  such  as  a piece  of  plaster  falling  from  the  ceiling,  is  the 
damage  apparent.  For  such  a failure,  where  the  plaster  actually 
breaks  loose  and  falls,  it  is  necessary  not  only  that  it  be  cracked  at 
its  surface,  but  also  that  its  hold  on  the  lath  structure  be  cracked. 

Patterns  of  plaster  failures  caused  by  sonic  booms  in  the  White 
Sands  tests  (30]  are  shown  in  Table  23.  Note  from  the  table  that 
almost  all  the  observed  damage  consisted  of  hairline  cracks  and  the 
extension  of  existing  cracks.  The  table  indicates  two  types  of  failure 
for  plaster,  diaphragm  and  racking.  In  diaphragm  failure  the  wall  or 
ceiling  bellies  out  due  to  the  sonic  boom,  bending  the  plaster  on  its 
surface.  In  racking  failure  the  adjacent  walls  and  ceiling  lean  forward 
slightly,  trying  to  push  the  wall  into  a parallelogram  shape.  This 
results  in  a tearing  action  on  the  plaster  at  the  corners  of  the  room 
and  near  doors  and  widows.  As  seen  from  the  table,  diaphragm 
failures  and  racking  failures  occur  at  about  the  same  overpressures. 
Since  racking  failures  are  much  harder  to  treat  analytically  and  data 
on  them  are  scarce,  we  treat  diaphragm  failures  only  here.  The 
results  should  also  extend  to  racking  failures,  however,,  on  the  basis 
of  the  similarity  in  failure  overpressures. 
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We  performed  calculations  of  the  probability  of  plaster  cracking 
from  sonic  boom  for  the  six  plaster  configurations  shown  in  Table  24. 
The  entries  in  the  table  for  Numbers  1 and  2 are  calculated  static 
pressures  fcr  the  ceiling  configuration  indicated.  The  first  ceiling 
has  low  strength  plaster  with  a tensile  strergth  of  100  psi,  and  the 
second  has  high  strength  plaster  with  a tensile  strength  of  350  psi. 

The  four  plaster  wall  configurations  shown  represent  experimental 
results  from  pressurizing  wall  sections  in  laboratories  [51,  52]. 
Plaster  configurations  3 and  4 failed  at  lower  overpressures  than 
would  be  expected  for  a typical  external  wall  in  a home,  since  th^se 
configurations  represented  party  walls  rather  than  outside  walls. 
Configurations  5 and  6 represent  the  most  susceptible  configurations 
of  external  walls,  those  of  wood  frame  houses.  The  fact  that  the 
failure  pressures  in  Table  24  are  so  low  shows  the  reason  why  plas- 
ter often  fails  in  ordinary  occurrences,  such  as  slamming  a door 
or  hammering  a nail  in  the  wall. 

A thorough  search  was  made  for  data  on  the  probability  density 
function  of  the  breaking  strength  of  plaster.  We  reviev/ed  the  litera- 
ture and  contacted  manufacturers,  trade  associations,  and  the 
National  Bureau  of  Standards.  Unfortunately  no  probaoility  density 
function  data  on  plaster  were  available,  We  were,  however,  able  to 
find  sufficient  data  o.i  a similar  material,  mortar  [37],  and  found 


Table  24.  Static  Failure  Pressures  of  Plaster  Configurations 
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its  distribution  to  be  lognormal.  The  fact  that  mortar  has  a lognormal 
distribution  allows  us  to  apply  the  response  probability  density  function 
technique,  as  we  did  for  glass,  ir  treating  brick-mortar  systems  later 
in  this  report.  On  the  basis  of  the  structural  similarity  of  plaster  and 
mortar,  we  represent  the  probability  density  function  of  plaster  here 
with  a lognormal  probability  density  function  having  the  variance  of  its 
logarithm  the  same  as  the  variance  of  the  logarithm  of  the  breaking 
pressure  of  mortar.  The  variance  of  the  logarithm  of  the  breaking 
pressure  of  mortar  was  found  to  be  . 0324;  this  figure  was  thus  also 
used  here  as  the  variance  cf  the  logarithm  of  the  breaking  pressure  of 
plaster. 

The  response  probability  density  function  technique  was  used  in 
analyzing  the  plaster  configurations  shown  in  Table  24.  The  resulting 
estimates  of  the  probability  of  breakage  as  a function  of  nominal 
overpressure  are  shown  in  Figure  30.  Note  that  the  breakage  proba- 
bilities indicated  are  comparable  to  those  of  glass. 

In  calculating  the  breakage  probabilities  of  the  ceiling  configur- 
ations 1 and  2,  different  values  of  the  mean  and  variance  of 
logio  (Pe^Pf)  were  used  than  were  used  for  the  walls.  This  is  done 
to  account  for  the  transmission  loss  of  the  boom  in  passing  thvough 
the  attic  structure.  For  ceilings  the  mean  of  logj0(pe/pf)  ' --1609 


Breakage  Probabilities  for  Plaster  for  All  Flight  Paths 
Equally  Likely 


and  the  variance  of  log^Q  (pe/pf)  = .0029.  For  the  wall  configurations 
the  mean  and  variance  of  logjQ  (pe/pf)  which  were  used  for  windows 
were  again  used.  The  dynamic  amplification  factor  distribution  found 
for  the  large  window,  Window  4,  was  used  with  the  plaster  calculations, 
since  this  window  is  comparable  in  size  to  a wall  or  ceiling.  On  the 
basis  of  the  above  procedures  the  graphs  snown  in  Figure  30  resulted. 

***#*#&**## 

In  summary,  the  estimated  breakage  probabilities  for  plaster, 
shown  in  Figure  30.  are  similar  to  the  results  found  for  glass.  As 
in  the  case  of  glass  there  is  a great  deal  of  variation  in  the  probab- 
bility  of  breakage  for  the  various  configurations. 
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vn.  BRIC-A-BRAC 

For  the  purpose  of  this  study  bric-a-brac  is  defined  as  an 
assortment  of  miscellaneous  ornamental  articles,  which  sit  on  sur- 
faces such  as  shelves  or  tables.  Typical  examples  would  be  figurines, 
ashtrays,  clocks,  cup  and  saucer  sets,  or  candlesticks.  Not  included 
in  this  definition  are  objects  which  hang  from  walls  such  as  pictures 
and  mirrors.  The  exclusion  of  hanging  objects  from  this  category  is 
in  keeping  with  earlier  sonic  boom  claims  categorization  work  [32  ] . 

By  its  very  nature  this  type  of  sonic  boom  damage  is  very  diffi- 
cult to  analyze.  The  objects  which  make  up  this  category  vary  by 
orders  of  magnitude  in  the  force  needed  to  overturn  or  slide  them. 
They  also  vary  by  orders  of  magnitude  in  their  probability  of  breaking 
once  they  fall  from  their  supporting  surface.  For  one  extreme, 
imagine  a heavy  steel  disc  paperweight  sitting  on  a rubber  mat  in 
the  middle  of  a table.  For  the  other  extreme,  imagine  a delicate 
saucer  whose  owner  has  balanced  it  on  edge,  precariously  leaning  on 
the  wall  of  a narrow  knick-knack  shelf,  with  a concrete  floor  below. 
The  steel  paperweight  would  be  immune  to  any  sonic  boom  ever 
created,  both  in  the  matter  of  falling  and  also  of  breaking  if  it  did 
fall.  The  precariously  balanced  saucer,  however,  would  probably 
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break  with  any  small  disturbance  such  as  a door  slamming,  a footfall, 
or  a wind  from  opening  a window. 

The  above  examples  point  out  three  key  factors: 

. There  is  no  "typical"  piece  of  bric  -a-brac. 

. Pieces  of  bric-a-brac  vary  tremendously  in  their 
strength. 

. The  susceptibility  of  pieces  of  bric-a-brac  to  overturning 
or  sliding  off  their  supports  depends  to  a large  extent  on 
how  their  owner  places  them. 

It  is  apparent  from  the  above  discussion  that  bric-a-brac  do  not 
lend  themselves  well  to  a stress  analysis,  such  as  was  used  for  glass 
windows.  The  determination  of  the  probability  of  breakage  of  bric-a- 
brac  will  have  to  be  based  instead  on  claims  data  and  on  observations 
of  bric-a-brac  breakage  during  severe  booms  at  White  Sands.  Our 
analysis  extrapolates  from  these  measured  breakage  probabilities  by 
assuming  a lognormal  effective  factor  of  safety  holds,  such  as  was 
observed  for  window  glass. 

The  two  measured  breakage  probabilities  for  bric-a-brac  are 
indicated  by  the  circles  on  Figure  31.  The  point  at  2 psf  is  from  the 
Edwards  Air  Force  Base  claims  data  [32  }.  The  point  at  38  psf  is 
from  the  White  Sands  experiment. 

The  break:  ze  probability  shown  at  2 psf  of  8.8  x 10"^  was 
calculated  by  observing  the  ratio  of  bric-a-brac  claims  to  window 
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Breakage  Probabilities  for  Bric-a-Brac  for 
All  Flight  Paths  Equally  Likely 


claims  from  the  data  in  Table  25  to  be  3/58.  We  also  estimated  that 
there  are  roughly  twice  the  number  of  bric-a-brac  items  in  an  average 
home  as  there  are  windows,  based  on  a count  in  the  author's  own  home. 
We  assumed  the  breakage  probability  of  a typical  window  population  to 
be  3.4  x 10 "5,  as  found  in  our  previous  studies.  Thus  multiplying 
these  gave  (3.4  x 10"^)  x (3/58)  x (1/2)  = 8.8  x 10-?  as  the  bric-a-brac 
breakage  probability. 

The  point  at  38  psf  indicating  a breakage  probability  of  0.2  per- 
cent was  derived  from  the  White  Sands  data  [30].  It  was  observed 
that  2 out  of  an  estimated  900  pieces  of  brie  -a-brac  in  the  test  houses 
were  broken  by  a sonic  boom  with  an  overpressure  of  38  psf. 

In  Figure  31  these  two  points  are  connected  by  a double  line  on 
a lognormal  probability  scale.  This  is  the  same  scale  that  was  used 
earlier  in  Figure  26  to  display  the  probability  of  breakage  of  glass  as 
a function  of  sonic  boom  nominal  overpressure.  The  graphs  for 
glass  breakage  are  shown  with  single  lines  on  the  same  graph  for 
comparison.  It  can  be  observed  from  the  very  steep  slope  of  the 
bric-a-brac  line  shown  that  the  standard  deviation  of  the  effective 
factor  of  safety  is  much  higher  for  bric-a-brac  than  it  is  for  windows. 
This  is  as  expected  since,  as  discussed  above,  bric-a-brac  has  a 
much  greater  structural  variation  than  does  glass.  The  bric-a-brac 
line  indicates  that  because  of  the  very  wide  variation  of  strength  and 
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support  condition  there  would  be  a breakage  probability,  comparable 
to  the  less  susceptible  glass,  even  at  very  low  nominal  overpressures, 
such  as  1 psf . By  the  same  token,  however,  it  is  unlikely  that  more 
than  about  1 percent  of  the  bric-a-brac  boomed  would  be  broken  even 
for  the  highest  overpressures. 


* * ******** 


In  summary,  bric-a-brac  breakage  is  generally  less  than  glass 
for  the  normal  range  of  sonic  boom  nominal  overpressures.  Estimates 
of  bric-a-brac  breakage  probabilities  are  given  by  the  line  in  Figure  31. 
However,  it  must  be  emphasized  that  the  breakage  probabilities  for 
bric-a-brac  exhibit  large  variations  depending  on  the  strength  of  the 
individual  items  and  their  support  conditions. 
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VIII.  BRICK 


Brick  is  a far  stronger  material  than  glass  or  plaster  and  is 
generally  not  susceptible  to  sonic  boom  damage.  For  the  sake  of 
completeness,  however,  we  are  also  considering  brick  in  this  report. 

Common  bricks  are  typically  manufactured  by  heating  rectangu- 
lar prisms  of  clay  or  shale  in  a kiln.  The  processes  used  include  the 
dry-press  process,  the  stiff-mud  process,  and  the  soft-mud  process 
[53] . In  the  dry-press  process  relatively  dry  clay  mixes  are  used 
and  the  molding  of  bricks  is  done  at  high  pressure.  This  results  in 
all  3ix  surfaces  being  smooth  and  even  with  good  uniformity  among 
the  bricks.  In  the  stiff-mud  process  the  clay  mixture  is  more  moist 
and  it  is  extruded  by  a machine  into  a long  ribbon  with  the  cross- 
section  of  a brick.  The  ribbon  is  then  cut  by  wires  into  brick-length 
pieces,  leaving  rough  surfaces  at  the  wire  cuts,  n the  soft -mud 
process  the  brick  is  formed  from  a wet  mix  of  clay  under  only  slight 
pressure. 

Various  forms  of  kilns  are  used  for  firing  the  brick.  In  tunnel 
kilns  the  bricks  are  placed  on  cars  and  move  slowly  through  the  tun- 
nel gradually  reaching  a temperature  of  about  2000°F,  and  then  they 
are  gradually  cooled  as  they  emerge.  The  tunnel  kiln  method  of 
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firing  results  in  very  uniform  bricks.  Other  \inds  of  kilns,  however, 
results  in  more  variation  among  bricks  depending  on  where  in  the 
kiln  the  brick  was  fired.  These  other  kilns  include  scove  kilns, 
round  kilns,  continuous  kilns,  and  permanent  kilns.  The  brick  may 
be  classified  as  "top  brick"  from  the  top  of  the  kiln,  "body  brick" 
from  the  center,  or  "bench  brick"  from  the  arches  of  the  kiln. 
Variations  in  strength  of  bricks  have  practically  no  effect  as  far  as 
sonic  boom  loading  is  concerned;  the  weak  point  in  the  brick  struc- 
ture is  almost  always  the  mortar. 


Mortar  consists  of  cement,  lime,  and  sand  combineu  in  vari- 
ous proportions.  The  proportions  used  result  in  mortars  with  desig- 
nations as  M,  S,  N,  O.  and  K as  shown  in  Table  26.  In  general  the 
compressive  strength  of  the  mortar  increases  with  increasing  pro- 
portions of  cement.  There  is  a trade-off,  however,  since  mortars 
with  more  lime  and  less  cement  tend  to  be  easiei  f »r  bricklayers  to 
apply. 

For  the  purpose  of  this  study  we  gathered  data  on  the  strength 
of  mortars  when  used  in  wah  assemblies  subjected  to  uniform  load- 
ing (37,54,55].  An  analysis  of  the  data,  covering  types  M,  N,  O, 
and  S mortars  showed  the  probability  density  function  of  the  strengths 
to  be  lognormaily  distributed  with  a mean  modulus  of  rupture  of  45.  6 
psi  and  a standard  deviation  of  16.  9 psi.  The  distribution  of  the 
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logarithms  of  the  strengths  had  a mean  of  1.  6252  and  a variance 
of  . 0324.  On  the  basis  of  these  data  probability  of  failure  of  brick 
structures  was  calculated  for  sonic  boom  loading  using  the  response 
probability  density  function  technique. 

The  first  calculation  presented  here  should  lay  to  rest  once 
and  for  all  any  fears  that  people  may  have  of  sonic  booms  knocking 
down  brick  walls  of  homes.  Table  27  from  the  nuclear  effects 
literature  [56]  shows  that  it  takes  a pressure  of  7 psi  (1008  psf)  to 
knock  down  a brick  wall  of  a house.  The  most  intense  sonic  booms 
deliberately  created  for  test  purposes  were  of  the  order  of  100  psf. 
Even  if  one  of  these  acts  on  the  brick  wall  the  breakage  probability 
is  1.3  x 10"15.  With  an  astronomical  probability  like  this  we  can 
say  that  for  all  practical  purposes  there  is  no  chance  of  knocking 
down  a brick  wall  of  a home. 

The  question  remains,  however,  of  the  probability  of  damage 
to  other  parts  of  more  susceptible  brick  structures;  fcr  instance, 
cracking  mortar  on  free-standing  brick  walls.  The  probability  of 
such  occurrences  will  be  examined  here  using  National  Bureau  of 
Standards  data  on  tests  of  brick  walls  [57]. 

The  four  brick  wall  configurations  analyzed  are  shown  in 
Table  28.  All  the  brick  walis  shown  were  tested  by  National  Bureau 
of  Standards  using  an  air  bag  technique  to  create  the  pressure  load. 
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Table  27.  Conditions  of  Failure  of  Peak 
Overpressure -Sensitive  Elements 


Structural  Element 

Foilwre 

Approximate 
Incident  Host 
Overpressure 
(psi) 

Gfose  windows,  largo  and 

Hood. 

Skattarinp  usuully, 
occasional  frame 
failure. 

0.5-I.0 

Corrugated  asbestos  siding 

Shattering. 

1. 0-2.0 

Corrugated  stool  or  olumi- 
inum  paneling. 

Connection  foilure 
followed  by  buclc- 
Bng. 

1. 0-2.0 

Wood  siding  panels,  stand- 
ard house  construction. 

Usually  failure  occurs 
at  the  noin  con* 
nections  allowing  a 
whole  panel  to  be 
blown  in. 

1. 0-2.0 

Concrete  or  cinder-block 
wall  panels  I in.  or  12 
in.  thick  (not  reinforced). 

Shotftring  of  the  wall. 

2.0-3.0 

•rtefc  WoU  panel,  8 in.  or 
12  in.  thick  (not  rein- 
forced). 

Shearing  and  flexure 
failures. 

2.0-8.0 

Table  28.  Static  Failure  Pressures  of  Brick  Wall  Configurations 


All  the  walls  were  constructed  of  bricks  with  a nominal  width  of 
4 inches.  The  static  failure  points  indicated  were  when  the  first 
cracking  of  the  mortar  was  observed.  In  each  case  two  specimens 
were  tested;  the  value  shown  in  Table  28  is  the  lower  of  the  two 
static  failure  pressures.  All  tests  were  performed  on  4-foot  wide 
by  8 -foot  high  test  panels.  The  mortar  referred  to  in  the  table  as 
1:1:4  consisted  (by  volume)  of  one  part  of  Type  1 portland  cement, 
one  part  hydrated  lime,  and  four  parts  sand.  The  mortar  referred 
to  in  the  table  as  "high-bond:t  contained  one  part  (1  cubic  foot)  port- 
land  cement,  one  part  (1  cubic  foot)  fine  limestone  (passing  a No.  200 
sieve),  four  parts  (4  cubic  feet)  masonry  sand,  and  four  gallons  of 
a polyvinyiidene  chloride  additive  called  Sarabond. 

The  breakage  probabilities  for  sonic  boom  loading  of  the  four 
brick  wall  configurations  were  calculated  using  the  response  proba- 
bility density  function  technique.  The  results  are  shown  in  Figure  32. 
Note  that  the  breakage  probabilities  shown  are  generally  lower  than 
those  for  glass  except  for  the  17.3  psf  wall,  which  was  not  made 
from  high -bond  mortar.  The  results  from  the  four  wall  cases  shown 
and  the  house  example  make  it  appear  extremely  unlikely  that  any 
damage  to  brick  structures  would  occur  during  supersoric  overflights. 
It  would  be  very  unusual  for  anyone  to  have  a free-standing  4-inch 
brick  wall,  and  only  one  made  from  lower  bond  mortar  would  be  more 
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IX.  CONCLUSIONS 


The  analysis  presented  here  has  generally  confirmed  the 
findings  from  sonic  boom  claims  experience.  The  most  likely 
elements  to  break  are  windows,  plaster,  and  bric-a-brac,  and  their 
likelihood  of  breaking  increases  with  increasing  overpressure. 

The  combined  graph  for  the  probability  of  breakage  of  glass, 
plaster,  bric-a-brac,  and  brick  is  shown  in  Figure  33.  It  presents 
the  results  that  had  been  shown  earlier  for  the  separate  materials, 
except  that  the  bric-a-brac  curve  has  been  broadened  to  indicate  a 
range  of  probabilities  at  each  overpressure  typical  of  other  mate- 
rials. The  curve  for  the  brick  wall  with  low  strength  mortar  was 
eliminated,  since  it  is  not  a very  representative  structure;  it  was 
merely  a sample  in  a laboratory  test  rather  than  a real  wall  in  use. 

For  sonic  booms  with  a nominal  overpressure  of  1 psf,  with 
all  flight  paths  equally  likely,  the  breakage  probability  ranges  for 
the  structural  elements  were  as  follows: 

. Windows:  4 x 10"®  to  3 x 10"® 

. Plaster:  3 x 10"^  to  5 x 10"^ 

Bric-a-Brac:  1 x 10"®  to  1 x 10_® 

. Brick  Walls:  **  x 10"^  to  1 x 10'^. 
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On  the  basis  of  the  results  of  the  study  various  measures  can 
be  suggested  for  reducing  the  probability  of  structural  damage  from 


supersonic  overflights.  These  measures  are  summarized  in  Table  29 
and  listed  below: 

1.  In  planning  flights — reducing  nominal  overpressures 

2.  In  designing  aircraft  — 

(1)  Increasing  altitude  capability 

(2)  Decreasing  weight 


3.  In  designing  window  installations  — 

(1)  Reducing  the  width  to  thickness  ratio 

(2)  Reducing  the  lite  size 

(3)  Replacing  cracked  panes 

(4)  Avoiding  scratches  during  installation  and  cleaning 

(5)  Placing  windows  so  that  they  face  away  from  the 
flight  path  of  supersonic  aircraft 

(6)  Mounting  windows  carefully  to  avoid  stress  con- 
centrations from  direct  contact  with  the  glazer's 
points  or  looseness  of  the  window  in  the  sash 

4.  In  designing  plaster  installations  — 

(1)  Stiffening  walls  and  ceilings 

(2)  Using  higher  strength  plaster 

(3)  Mixing  plaster  thoroughly  with  correct  proportion 
of  water  during  installation 

5.  In  placement  of  bric-a-brac — avoiding  precarious 

positioning  of  items  near  edges 
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6. 


In  designing  brick  walls  — 

(1)  Using  thicker  walls 

(2)  Laying  of  bricks  with  high-bond  mortar. 


I 


i 


'€ 


In  conclusion,  this  report  has  estimated  the  probabilities  of 
structural  damage  of  various  susceptible  elements  using  the  response 
probability  density  function  technique.  These  probabilities  were 
found  to  vary  widely  with  the  specific  material  configuration,  but 
to  consistently  increase  with  increasing  nominal  overpressures. 

For  all  materials  in  good  condition,  breakage  probabilities  were 
found  to  be  below  . 001  for  a nominal  overpressure  of  1 psf.  By 
using  appropriate  measures  in  aircraft  design,  flight  planning,  and 
material  installation  the  probability  of  sonic  boom  structural 
damage  can  be  significantly  reduced. 
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